Final Analysis of Proton Form Factor Ratio Data at = 4.0, 4.8 and 5.6 GeV^ 
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polarization transfer method at Jefferson Lab have revolutionized the understanding of nucleon 
structure by revealing the strong decrease of R with momentum transfer for > 1 GeV^, in 
strong disagreement with previous extractions of R from cross section measurements. In particular, 
the polarization transfer results have exposed the limits of applicability of the one-photon-exchange 
approximation and highlighted the role of quark orbital angular momentum in the nucleon structure. 
The GEp-II experiment in Jefferson Lab's Hall A measured R at four values in the range 3.5 
GeV^ < < 5.6 GeV^. A possible discrepancy between the originally published GEp-II results 
and more recent measurements at higher motivated a new analysis of the GEp-II data. This 
article presents the final results of the GEp-II experiment, including details of the new analysis, 
an expanded description of the apparatus and an overview of theoretical progress since the original 
publication. The key result of the final analysis is a systematic increase in the results for R, 
improving the consistency of the polarization transfer data in the high-Q^ region. This increase is 
the result of an improved selection of elastic events which largely removes the systematic effect of 
the inelastic contamination, underestimated by the original analysis. 
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I. INTRODUCTION 

The electromagnetic form factors (FFs) of the nucleon 
have been revived as a subject of high interest in hadronic 
physics since a series of precise recoil polarization mea- 
surements of the ratio of the proton's electric (G^) and 
mag netic {Glj) FFs [J, i in Jefferson Lab's Hall A es- 
tablished the rapid decrease with momentum transfer 
of R = iipG^^/G^j, where /ip is the proton's magnetic 
moment, for 0.5 GeV^ < < 5.6 GeV^. These mea- 
surements disagreed strongly with previous extractions 
of G^ from cross section data [3] using the Rosenbluth 
method [i*], which found /ipG^/G^^ ~ 1. Subsequent in- 
vestigations of both experimental techniques, including a 
novel "Super- Rosenbluth" measurement using ^H{e,p)e' 
cross section measurements to reduce systematic uncer- 
tainties [1| , found no neglected sources of error in either 
data set, pointing to incompletely understood physics as 
the source of the discrepancy. 

Theoretical investigations of the discrepancy have fo- 
cused on higher-order QED corrections to the cross sec- 
tion and polarization observables in elastic ep scattering 
[^, 0], including radiative corrections and two-photon- 
exchange (TPEX) effects. The amplitude for elastic 
electron-proton scattering involving the exchange of two 
or more harc0 photons cannot presently be calculated 
model-independently. In the region of the discrep- 
ancy, model calculations of TPEX [8, 9] predict relative 
corrections to both the cross section and polarization 
observables that are typically at the few-percent level. 
At large Q^, the sensitivity of the Born (one-photon- 
exchange) cross section to G^ becomes similar to or 
smaller than the sensitivity of the measured cross sec- 
tion to poorly-known TPEX corrections, obscuring the 
extraction of G^. On the other hand, the polarization 
transfer ratio R defined in equations ([l} is directly pro- 
portional to G^/G^^, such that the extraction of G^ is 
much less sensitive to corrections beyond the Born ap- 
proximation. For this reason, a general consensus has 
emerged that the polarization transfer data most reliably 
determine G|, at large Q^. Nonetheless, active exper- 
imental and theoretical investigation of the discrepancy 
and the role of TPEX continues 1^ . Owing to the lack of 
a model-independent theoretical prescription for TPEX 
corrections, precise measurements of elastic ep scattering 
observables sensitive to TPEX effects continue to play an 
important role in the resolution of the discrepancy. 

The revised experimental understanding of the proton 
form factors led to an onslaught of theoretical work. The 
constancy of the Rosenbluth data for G^/G^^ was con- 
sistent with a "p recocious" onset of pQCD dimensional 
scaling laws valid for asymptotically large Q^, an 
interpretation which had to be abandoned in light of 



the polarization data. The decrease of R with was 
later interpreted in a pQCD-scaling framework including 
higher-twist corrections [l^l, demonstrating the impor- 
tance of quark orbital angular momentum in the inter- 
pretation of nucleon structure. The relations between 
nucleon form factors and Generalized Parton Distribu- 
tions (GPDs) have placed this connection on a more 
quantitative footing [isl - fTsj . Furthermore, the GPD- 
form factor sum rules have been used to derive model- 
independent representations of the nucleon transverse 
charge and magnetization densities as two-dimensional 
Fourier transforms of the Dirac (Fi) and Pauli (F2) form 
factors [1^. In the context of calculations based on 
QCD's Dyson-Schwinger Equations (DSEs) [13, El, the 
form factor data are instrumental in elucidating the dy- 
namical interplay between the nucleon's dressed-quark 
core, diquark correlations, and the pseudoscalar meson 
cloud [19]. Recent measurements of the neutron form 
factors at large 20, 21] have enabled for the first 
time a detailed flavor-decomposition [22] of the form fac- 
tor data, leading to new insights. In addition, the form 
factor data have been interpreted within a large number 
of phenomenological models; a recent review of the large 
body of theoretical work relevant to the nucleon FFs is 
given in 1^, and a current overview is given in section 
HVBl of this work. 

The recoil polarization method exploits the relation 
between the transferred polarization in elastic ep scatter- 
ing and the ratio G^/G^^. In the one-photon-exchange 
approximation, the polarization transferred to recoiling 
protons in the elastic scattering of longitudinally po- 
larized electrons by unpolarized protons has longitudi- 
nal {Pg) and transverse (Pt) components in the reaction 
plane given by ^ [13] 



Pt = -hP, 
Pt = hP, 



2e(l - e) r 
T l + ^r2 



1 + 



Lr2 



(1) 



Pt_/r(l + e) _ R 
P, 



2e 



where h — ±.1 \s the electron beam helicity, P^. is the 
beam polarization, r = Q"^ / AM^, Mp is the proton mass, 

and e = [l + 2{1 + t) ta.n^ {ee/2)]~^ , with Be the elec- 
tron scattering angle in the proton rest (lab) frame, cor- 
responds to the longitudinal polarization of the virtual 
photon in the one-photon-exchange approximation. 

Recent measurements from Jefferson Lab's Hall C 
[2^ extended the reach of the polarization transfer 
method to 8.5 GeV^. The published data from Hall A 
are well described by a linear dependence [1], 



i?= 1.0587 - 0.14265Q2 



(2) 



^ "Hard" in this context means that both exchanged photons carry 
an appreciable fraction of the total momentum transfer. 



with 0^ in GeV2, valid for > 0.4 GeV^. On the other 
hand, all three of the recent Hall C data points are at 
least 1.5 standard deviations above this line, including 
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the measurement at overlapping = 5.17 GeV^, which 
Ues 1.8a above equation Due to the strong, incom- 
pletely understood discrepancy between the Rosenbluth 
and polarization transfer methods of extracting G^/G^^ 
and the fact that the new Hall C measurements arc the 
first to check the reproducibility of the Hall A data us- 
ing a completely different apparatus in the region 
where the discrepancy is strongest, understanding pos- 
sible systematic differences between the experiments is 
important. 

This article reports an updated, final data analysis of 
the three higher- measurements of R from Hall A, orig- 
inally published in ^2] , along with expanded details of the 
experiment. To avoid confusion, a naming convention is 
adopted throughout the remainder of this article for the 
most frequently cited experiments: GEp-I for Ref. 
GEp-H for Ref. 0, the subject of this article, GEp-HI for 
Ref. \2§\ and GEp-27 for Ref. SectionlUpresents the 
kinematics of the measurements, an expanded descrip- 
tion of the experimental apparatus, and a comparison of 
the GEp-H and GEp-HI/GEp-27 experiments. Section 
IIHI presents the data analysis method, including the se- 
lection of elastic events, the extraction of polarization ob- 
servables, and the estimation and subtraction of the non- 
elastic background contribution. Section ITVl presents the 
final results of the experiment and discusses the impact 
of the revised data on the world database of proton elec- 
tromagnetic form factor measurements, in the context 
of the considerable advances in theory since the original 
publication. The conclusions and summary are given in 
section |Vl 



II. EXPERIMENT SETUP 

Table U shows the central kinematics of the measure- 
ments from the GEp-H experiment. The kinematic vari- 
ables given in Table H] are the beam energy , the scat- 
tered electron energy E'^,, the electron scattering angle 
9e, the scattered proton momentum pp, and the proton 
scattering angle 9p. 



The Continuous Electron Beam Accelerator at the 
Thomas Jefferson National Accelerator Facility (JLab) 
delivers a high quality, longitudinally polarized electron 
beam with ^100% duty factor. The beam energy was 
measured using the Arc and ep methods. The ep method 
determines the energy by measuring the opening angle 
between the scattered electron and the recoil proton in 
ep elastic scattering, while the Arc method uses the stan- 
dard technique of measuring a bend angle in a series of 
dipole magnets. The combined absolute accuracy of both 
methods is AE/E ~ 10~*, while the beam energy spread 
is at the 10~^ level. The nominal beam energy in this ex- 
periment was 4.6 GeV. The beam polarization was mea- 
sured by Compton and Moller polarimeters. Details of 
the standard Hall A equipment can be found in [27| and 
references therein. 

The hydrogen target cell used in this experiment was 
15 cm long along the beam direction. The target was 
operated at a constant temperature of 19 K and pressure 
of 25 psi, resulting in a density of about 0.072 g/cm^. To 
minimize the target density fluctuations due to localized 
heat deposition by the intense electron beam, a fast raster 
system consisting of a pair of dipole magnets was used to 
increase the transverse size of the beam in the horizontal 
and vertical directions. The raster shape was square or 
circular in the plane transverse to the beam axis. In this 
experiment, the raster size was approximately 4x4 mm^. 

Recoil protons were detected in the high resolution 
spectrometer located on the beam left (HRSL) ^2l\ . The 
HRSL has a central bend angle of 45°, and subtends a 
6.5 msr solid-angle for charged particles with momenta 
up to 4 GeV with ±5% momentum acceptance. Two ver- 
tical drift chambers measure the particle's position and 
trajectory at the focal plane. With knowledge of the 
optics of the HRSL magnets and precise beam position 
monitoring, the proton scattering angles, momentum and 
vertex coordinates were reconstructed with FWHM reso- 
lutions of ~2.6 (4.0) mrad for the in-plane (out-of-plane) 
angle, Ap/p ^ 2.6 x 10^"* for the momentum, and ^3.1 
mm for the vertex coordinate in the plane transverse to 
the HRSL optical axis. 



A. Experimental Apparatus 

The GEp-II experiment ran in Hall A at Jefferson Lab 
during November and December of 2000. A polarized 
electron beam was scattered off a liquid hydrogen tar- 
get. Hall A is equipped with two High Resolution Spec- 
trometers (HRS) [27}, which are identical in design. In 
this experiment, the left HRS (HRSL) was used to detect 
the recoil proton, while the right HRS (HRSR) was used 
to detect the scattered electron at the lowest of 3.5 
GeV^. For the three highest points at 4.0, 4.8 and 5.6 
GeV^, electrons were detected by a lead-glass calorime- 
ter. The focal plane of the HRSL was equipped with a 
focal plane polarimeter (FPP) to measure the polariza- 
tion of the recoil proton. 



1. Focal Plane Polarimeter 

The central instrument of this experiment was the Fo- 
cal Plane Polarimeter (FPP) Q, installed in the focal 
plane of HRSL. The FPP measures the transverse polar- 
ization of the recoil proton. The protons are scattered in 
the focal plane region by an analyzer, as shown in Fig- 
ure [TJ If the protons are polarized transverse to their 
momentum direction, an azimuthal asymmetry results 
from the spin-orbit interaction with the analyzing nuclei. 

The FPP has been described in detail in [l[, so only 
a brief summary of its characteristics will be given here. 
The only significant difference in the configuration of the 
FPP between the GEp-I and GEp-II experiments was a 
change of the analyzer material from carbon to polyethy- 
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TABLE I. Central kinematics of the GEp-II experiment. The central value is defined by the central momentum of the left 
HRS in which the proton was detected, e is the parameter appearing in equation ([T]), i5e is the beam energy, E'^ is the scattered 
electron energy, 9^ is the electron scattering angle, pp is the proton momentum, Op is the proton scattering angle, x is the 
central spin precession angle, Pe is the beam polarization, and Real is the distance from the target to the calorimeter surface. 
At the lowest of 3.5 GeV^, the scattered electron was detected in the right HRS. 



Nominal (GeV^) 


e 


Ee (GeV) 


£e (GeV) 


de n 


Pp (GeV) 


Op n 


xil 


Pe (%) 


Real (m) 


3.5 


0.77 


4.61 


2.74 


30.6 


2.64 


31.8 


241 


70 


HRSR 


4.0 


0.71 


4.61 


2.47 


34.5 


2.92 


28.6 


264 


70 


17.0 


4.8 


0.59 


4.59 


2.04 


42.1 


3.36 


23.8 


301 


73 


12.5 


5.6 


0.45 


4.60 


1.61 


51.4 


3.81 


19.3 


337 


71 


9.0 



Rear straw chambers 



Front straw chambers 

,/\ 




CH2 analyzer 

FIG. 1. (color online) Layout of the Focal Plane Polarimeter. 



lene. During GEp-I, the analyzer consisted of four doors 
of carbon which could be combined to produce a maxi- 
mum thickness of 51 cm. For cost, safety and efficiency 
reasons, carbon is ideal for measuring proton polarization 
with a momentum up to 2.4 GeV, which was sufficient 
for GEp-I. For GEp-II, the maximum proton momentum 
was 3.8 GeV. At this momentum, the analyzing power of 
carbon, which contributes to the size of the asymmetry, 
and therefore to the size of the error bar, decreases signif- 
icantly. An experiment was carried out at the Laboratory 
for High Energy (LHE) at the Joint Institute for Nuclear 
Research (JINR) in Dubna, Russia to find an optimal an- 
alyzing material and its thickness for protons at 3.8 GeV 
[28| . Polyethylene, a compound of carbon and hydrogen, 
was found to increase the analyzing power relative to car- 
bon as shown in Figured] A stack of eighty 2.5 cm-thick 
plates, each 58 cm deep along the direction of incident 
protons, was installed between the unused, opened, doors 
of the carbon analyzer, as shown in Figure[31 This 58 cm 
thickness was used for the — 3.5 GeV^ kinematics. 
For the three higher kinematics, an additional stack 
of polyethylene with a thickness of 42 cm was installed 
on a rail just upstream of the 58 cm stack to give a total 
thickness of 100 cm. 
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FIG. 2. (color online) Analyzing power versus scattering angle 
in the analyzer for graphite (49.5 cm-thick) and polyethylene 
(58 cm-thick) at a proton momentum of 2.64 GeV, corre- 
sponding to — 3.5 GeV^. 




FIG. 3. (color online) Stack of polyethylene plates for the 
analyzer. The dimensions shown on the plate are for the 
58 cm (42 cm) stack, and were chosen to match the envelope 
of elastically scattered protons in HRSL. 



2. Electron detection at Q = 3.5 GeV 

For the measurement at = 3.5 GeV^, the electron 
was detected in the high resolution spectrometer located 
on the beam right (HRSR). The trigger was defined by a 
coincidence between an electron in HRSR and a proton in 
HRSL. The precise measurement of the scattered electron 
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kinematics using a high-resolution magnetic spectrome- 
ter provides for an extremely clean selection of elastic ep 
events with cuts on the reconstructed mis sing energy and 
momentum, as shown in Fig. 4.8 of Ref. [29| . 



3. Electron detection at > 4.0 GeV^ 



For the measurements at > 4.0 GeV^, a lead- 
glass calorimeter was used to detect electrons due to the 
larger electron solid angle compared to the proton solid 
angle defined by HRSL. The lead- glass blocks from the 
standard HRSR calorimeter were used to assemble this 
calorimeter along with some additional spare blocks. Fig- 
ure 0] shows a front and a side view of the calorimeter on 
its platform. The blocks of lead-glass, of cross-sectional 
area 15 x 15 cm^, were individually wrapped in one foil 
of aluminized mylar, and one foil of black paper, to avoid 
light leaks. Each block was then tested, and the current 
drawn in the phototube due to noise was found to be less 
than 100 nA for all blocks. The blocks were assembled 
in a rectangular array of 9 columns and 17 rows, requir- 
ing a total of 153 blocks. Most of the blocks, in green 
in Figure m were 35 cm long, corresponding to 13.7 ra- 
diation lengths. 37 blocks positioned on the edges of 
the calorimeter were only 30 cm long, corresponding to 
11.8 radiation lengths (in blue in Figure 2]). The highest 
electron energy, at — 4.0 GeV^, was 2.5 GeV. At this 
energy, the shower stops after 7.7 radiation lengths, so 
that the entire shower is contained in the block. Since 
only 147 blocks were available, 6 blocks were missing to 
form a complete rectangle. These were replaced by wood 
placeholder blocks at the corners of the detector (in red in 
FigureH]) . The active area of the calorimeter was 3.31 m^ . 

The blocks were placed in the steel support frame (1), 
and held together using wood plates (2). The front of 
the support was covered with a one inch thick aluminum 
plate (3) to absorb very low energy particles. The ensem- 
ble was lifted by the top steel plate (4) onto the platform 
(5) using the Hall A crane. Balance on the platform was 
maintained by the steel support legs (6). The ensemble 
was placed on wheels and moved with the help of the 
Hall A crane attached to the steel lifting frame (7). The 
calorimeter was placed at the distance from the target 
required to match the electron solid angle to that of the 
proton at each kinematic setting. The acceptance match- 
ing was only approximate, due to the complicated shape 
of the spectrometer acceptance. Overall, about 5% of 
elastic events with a proton detected in HRSL were lost 
due to acceptance mismatching. The Cherenkov light 
emitted by primary electrons and shower secondaries in 
the lead-glass was collected by Photonis XP2050 photo- 
multipliers (8), and the signals were digitized by LeCroy 
1881 integrating ADCs and LeCroy 1877 TDCs. The 
trigger for the measurements at > 4.0 GeV^ was de- 
fined by a single proton in the HRS, signaled by a coin- 
cidence of two planes of scintillators in the focal plane. 
For each single proton event in the left HRS, the ADC 



(2): 







35 cm block 



Front view 



Side view 



FIG. 4. (color) Design of the calorimeter used to detect the 
scattered electron. In the front view, the 2.54 cm-thick alu- 
minum plate in front of the blocks is not shown. See text for 
details. 



and TDC information from the calorimeter was read out 
for all blocks, and elastic events were selected offline by 
applying software cuts to the calorimeter data. 



B. Comparison to Hall C Experiments 

The GEp-II experiment shares several important fea- 
tures with GEp-III. Both experiments used magnetic 
spectrometers instrumented with Focal Plane Polarime- 
ters (FPPs) to detect protons and measure their polariza- 
tion, and large acceptance electromagnetic calorimeters 
to detect electrons in coincidence. The use of calorime- 
ters in both experiments was driven by the requirement of 
acceptance matching; at large and 9e, the Jacobian of 
the reaction magnifies the electron solid angle compared 
to the proton solid angle fixed by the spectrometer ac- 
ceptance. The drawbacks of this choice compared to elec- 
tron detection using a magnetic spectrometer are twofold. 
First, the energy resolution of lead-glass calorimeters is 
relatively poor, so that elastic and inelastic reactions are 
not well separated in reconstructed electron energy. Sec- 
ond, the signals in lead-glass from electrons and pho- 
tons of similar energies are indistinguishable, leaving one 
vulnerable to photon backgrounds from the decay of tt*^, 
which played an important role in the analysis of both 
experiments. 

The high-Q^ measurements of the GEp-III experiment 
[25} were carried out consecutively with the GEp-27 
experiment, a series of precise measurements of R at 
— 2.5 GeV^ [11] designed to search for effects be- 
yond the Born approximation, thought to explain the 
disagreement between Rosenbluth and polarization data 
. Using the same apparatus and analysis procedure as 
GEp-III, the results of GEp-27 [26] arc in excellent agree- 
ment with the GEp-I data from Hall A [1] at nearly iden- 
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tical Q^: as shown in Figure [T2l The background correc- 
tions to the GEp-27 data were neghgible after applying 
the cuts described in [H, [ll] . Similarly, electrons were 
detected in the HRSL in the GEp-I experiment, so that 
the selection of elastic events was practically background- 
free In the absence of major background corrections, 
the agreement between precise measurements at the same 
using different polarimeters and magnetic spectrome- 
ters limits the size of any potentially neglected systematic 
errors arising from sources other than background. 



The liquid hydrogen targets used in Halls A and C 
had radiation lengths of ~2%, leading to a significant 
Bremsstrahlung fiux across the target length, in addi- 
tion to the virtual photon flux due to the presence of the 
electron beam. The kinematics of tt*' photoproduction 
(7 -l-p — ?> tt'^ +p) near end point {E^ Ee) are very sim- 
ilar to elastic ep scattering at high energies {E^ ^ m,^), 
such that protons from 7 -I- p — ^ tt*^ -f p overlap with elas- 
tically scattered protons within experimental resolution. 
In the lab frame, asymmetric tt" decays with one photon 
emitted at a forward angle relative to the tt" momen- 
tum, carrying most of the energy, are detected with a 
high probability. At high energies and momentum trans- 
fers, the 7r° photoproduction cross section is observed to 
scale as for fixed Qcm 0, where s is the center- 
of-momentum (CM) energy squared and Qcm is the CM 
7r° production angle. In addition, the CM angular distri- 
bution is peaked at forward and backward angles. The 
goal of the GEp-III experiment was to measure to the 
highest possible Q^, given the maximum available beam 
energy of 5.71 GeV. At = 8.5 GeV^, the relatively 
high Q'^/s ratio, with Qcm e 129-143°, led to a TT°p:ep 
ratio of ~40:1. The severity of the background condi- 
tions required maximal exploitation of elastic kinematics 
to suppress the 7r° background. Even after all cuts de- 
scribed in [25[ , the remaining background was estimated 
at ~6% of accepted events. Given the large difference be- 
tween the signal and background polarizations, this level 
of contamination required a substantial positive correc- 
tion to R. 



In light of the improved understanding of the impor- 
tance of the 7r° background gained during the analysis 
of the GEp-III data, an underestimation of its effect in 
the GEp-II analysis was considered as a potential source 
of disagreement between the two experiments. There- 
fore, the GEp-II data for = 4.0, 4.8, and 5.6 GeV^ 
were reanalyzed to investigate the systematics of the tt" 
background. The data from GEp-II at = 3.5 QeV^ 
were not reanalyzed, since electrons were detected in the 
HRSR and the 7r° background was absent. The systemat- 
ics of this configuration were thus irrelevant to the com- 
parison between GEp-II and GEp-III at higher Q^. 



III. DATA ANALYSIS 



Elastic Event Selection 



Figure [5] shows a representative example of the proce- 
dure for isolating elastic events in the GEp-II data, at 
Q'^ — 4.8 GeV^. As described in [2] and |2g|, cuts were 
applied to the difference between the HRS and calorime- 
ter time signals (±4 ns at = 4.0 and 4.8 GeV^, 
and ±5 ns at — 5.6 GeV^) and the missing energy 

( E 



Ee+M^ 



~ y Pp ~^ ^p ~ ^calo 

to suppress random coincidences and low-energy inelastic 
backgrounds, respectivelj0. The remaining backgrounds 
from ^11(7, TT^p) and quasielastic A\{e,e'p) reactions 
in the target cell windows were rejected using the 
kinematic correlations between the electron and proton 
arms. The measured proton kinematics were used to 
predict the scattered electron's trajectory assuming 
elastic scattering, and then the predicted electron trajec- 
tory, defined by the polar scattering angle 6*6^-* and the 
azimuthal scattering angle ^l^-* (where ^^'^ denotes the 
value predicted from the measured proton kinematics), 
was projected from the measured interaction verte>|f| to 
the surface of the calorimeter. 

Figures [SJa) and [IJb) show the horizontal (Ax) and 
vertical (Ay) differences between the measured shower 
coordinates at the calorimeter and the coordinates cal- 
culated from the measured proton kinematics assum- 
ing elastic scattering. Figure [5Kc) shows the difference 
Sp = Pp(Op) — Pp between the measured Pp and the mo- 
mentum required by elastic kinematics at the measured 
dp, given by 



Er. 



< 1000 MeV 



Pp{Op) 



2MpEe{Mp + Ee) COS Op 

2MpEe + E'i sin^ 6'„ 



(3) 



In each panel of Figure [5l the distribution of the plotted 
variable is shown before and after applying cuts (illus- 
trated by vertical lines) to both of the other two variables, 
which most nearly corresponds to the GEp-III analysis. 
In addition, the Ax (Ay) distribution is shown after ap- 
plying the cut to Ay (Ax), regardless of Sp, which most 
nearly corresponds to the selection of the original GEp-II 
analysis, in which no cut was applied to 6p. Each spec- 
trum exhibits a clear elastic peak near zero on top of 
a smooth background distribution. The background in 
the Ax and Ay spectra is dominated by 7r° photopro- 
duction events. The estimated background curves shown 
in panels (a) and (b) of Figure [S] were obtained using 
the polynomial sideband fitting method described in sec- 
tion |lirc2] The Sp cut clearly has significant additional 



^ The loose missing energy cut reflects the relatively poor energy 
resolution of lead-glass. 

^ The interaction vertex is defined as the intersection of the beam- 
line with the projection of the reconstructed proton trajectory 
on the horizontal plane. 
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FIG. 5. (color online) Elastic event selection at = 4.8 GeV^. The effects of cuts are shown for the horizontal calorimeter 
coordinate difference Aa; in panel (a), the vertical difference Ay in panel (b), and the proton momentum difference &p = 
Pp{Op) — Pp in panel (c). Solid vertical lines indicate the cut applied. Empty triangles show the distribution of events before 
applying any cuts. Filled squares show events passing the cuts on both of the other two variables. Empty circles in panels (a) 
and (b) show the Aa; (Ay) distribution of events passing the Ay (Ax) cut, regardless of 5p. In panels (a) and (b), the dashed 
and solid curves show the estimated background before and after the Sp cut. Dot-dashed vertical lines indicate the range of 
the elastic peak excluded from the fit to the background. 



background suppression power relative to Aa; and Ay 
cuts alone. In the 6p spectrum, the background distribu- 
tion is highly asymmetric about the peak, reflecting the 
fact that elastically scattered protons carry the highest 
kincmatically allowed momenta at a given Op. 

Since the two-body reaction kinematics are overdeter- 
mined, the method used to calculate Ax and Ay is not 
unique. In combination with the precisely known beam 
energy, the expected electron polar scattering angle o'f^ 
can be calculated from either the measured proton mo- 
mentum Pp , the measured proton scattering angle Op , or a 
combination of both. Different methods were used by the 
GEp-II and GEp-III data analyses to calculate Aa; and 
Ay. In the original GEp-II analysis, the calculation was 
formulated in terms of Cartesian components of the out- 
going particle momenta rather than polar and azimuthal 

(v) 

scattering angles. The effective 9e in the GEp-II ap- 
proach depends on both 9p and Pp. The exact equa- 
tions used can be found in Appendix D of Ref. [1^ . In 
the GEp-III analysis, 0^^' was calculated from pp, as de- 
scribed in [2^ . Both methods were tested in the present 
reanalysis. The Aa; and Ay distributions in Figures [SJa) 
and Eljb) were calculated using the GEp-II method, in 
order to demonstrate the background suppression power 
of the added Sp cut of Fig. [SJc) relative to the original 
analysis. For events selected by this cut, pp « Pp{Op), 
such that the Aa; values obtained from the GEp-II and 
GEp-III methods are equal up to detector resolution. 

A key difference between the GEp-III and GEp-II ex- 
periments is the dominant source of resolution in the 
variables used to select elastic events. The cell size of 
the GEp-II calorimeter was 15 x 15 cm^, compared to 
the 4x4 cm^ cell size of the GEp-III calorimeter. In 
GEp-II, the resolution of Aa; and Ay is dominated by 



the calorimeter coordinate measurement, and is therefore 
largely insensitive to the choice of proton variables used 
to calculate the expected electron angles. In GEp-III, on 
the other hand, the scattered electron angles were mea- 
sured with excellent resolution by the highly-segmented 
BigCal, such that the proton arm resolution was domi- 
nant. Given the kinematics of GEp-III and the angular 
and momentum resolution of the High Momentum Spec- 
trometer (HMS) in Hall C ^l|, the best Ax resolution 

was obtained by using pp to calculate Oi^^ . In the GEp-II 
analysis, the main practical difference between the two 
methods is the resulting background shape. In kinemat- 
ics for which the reaction Jacobian necessitates the use of 
a calorimeter for electron detection, the GEp-III method 
generally results in a wider and more asymmetric Aa; dis- 
tribution of the background, with inelastic events assum- 
ing predominantly negative Aa; values. In the GEp-III 
analysis, using o'f\pp) provided the best possible Aa; res- 
olution and a wider Aa; distribution of the background. 
In the GEp-II case, calculating Aa; using the GEp-III 
method spreads out the background without affecting the 
width of the elastic peak, thus reducing the background 
in the Ax spectrum with no 6p cut. After applying the Sp 
cut, however, the Ax distributions obtained from the two 
calculations are practically identical, and the choice be- 
comes arbitrary. As discussed in section lTlI Dl the results 
for R obtained with the Sp cut included do not depend 
on the method used to calculate Ax. The final reanalysis 
results were obtained with Ax and Ay calculated using 
the GEp-II method. 

The original GEp-II analysis applied a two-dimensional 
polygon cut to the correlated Ay versus Ax distribution. 
Using identical cuts to the original analysis, the pub- 
lished results 1^ were successfully reproduced. In the fi- 



8 



nal analysis, however, one-dimensional (rectangular) cuts 
were applied to Aa; and Ay, which simplifies the back- 
ground estimation procedure. For all three points, a 
cut of ±12(±16) cm was applied to Ax{Ay), centered at 
the midpoint between half-maxima on either side of the 
elastic peak, as in Figure \E\ The width of the cuts was 
chosen to be similar to the effective width of the poly- 
gon cut applied by the original analysis, and reflects the 
dominant contribution of the calorimeter cell size to the 
resolution of Ax and Ay. In addition, a cut of ±15 MeV, 
also centered at the midpoint between half-maxima of the 
elastic peak, was applied to 6p, as in Figure [SJc) . The 
width of the Sp cut was chosen to be ±3(t, where tr « 5 
MeV is the Sp resolution, which was roughly independent 
of the proton momentum in this experiment. While the 
difference in the selection of events from using a different 
shape of the Ax and Ay cuts is small, the Sp cut removes 
a rather substantial 6.0%, 7.3%, and 10.7% of events rel- 
ative to the original analysis for — 4.0, 4.8 and 5.6 
GeV^, respectively. 

While a fraction of the events rejected by the Sp cut 
are elastic, including events in the ep radiative tail and 
elastic events with Sp smeared by non-Gaussian tails of 
the HRS resolution, most of the rejected events are part 
of the background, and contribute very little to the sta- 
tistical precision of the data. Moreover, even real elastic 
events reconstructed outside the peak region of Sp do 
not meaningfully contribute to the accurate determina- 
tion of the form factor ratio, because such events are 
either (a) part of the radiative tail and therefore subject 
to radiative corrections that are in principle calculable Q 
but practically difhcult due to large backgrounds in the 
radiative tail region, or (b) have unreliable angle or mo- 
mentum reconstruction, which distorts the spin transport 
matrix of the HRS (see Ref. [S^j and section lTlI B 2l below) 
in an uncontrolled fashion. Therefore, the application of 
the Sp cut has benefits beyond mere background suppres- 
sion, as it also suppresses radiative corrections and the 
(potential) systematic effects of large angle or momentum 
reconstruction errors. The estimation of the background 
contamination and the background-related corrections to 
the polarization transfer observables are discussed in sec- 
tion IIII CI The next section discusses the procedure for 
the extraction of polarization observables from the "raw" 
asymmetries measured by the FPP. 



B. Extraction of Polarization Observables 

As detailed in [l|, [2^ , useful scattering events in the 
FPP were selected by requiring a good reconstructed 
track in both the front and rear straw chambers and 
requiring the scattering vertex Zdose, defined by the 
point of closest approach between incident and scattered 
tracks, to lie within the physical extent of the CH2 ana- 
lyzer. Events with polar scattering angles ?? < 0.5° were 
rejected, since at small angles comparable to the angu- 
lar resolution of the FPP, the azimuthal angle resolution 



diverges. Moreover, the small-angle region is dominated 
by multiple Coulomb scattering, which has zero analyz- 
ing power. 

1. Focal Plane Asymmetry 

Spin-orbit coupling causes a left-right asymmetry in 
the angular distribution of protons scattered by carbon 
and hydrogen nuclei in the CH2 analyzer of the FPP 
with respect to the transverse polarization of the incident 
protorQ. The measured angular distribution for incident 
protons with momentum p and transverse polarization 
components P^^^ and Py^^ for a beam helicity of ±1 
can be expressed a^ 

7V±(p, t?, ^) = Nf-^ [1 + {±AyPr'' + ci) cos ^+ 

(TAj,ff^^ + si)sin(^+ 

C2 cos(2v?) + 32 sin(2(p) + ...], (4) 

where is the total number of incident protons for 
beam helicity ±1, e{p,'d) is the polarimeter efficiency 
defined as the fraction of protons of momentum p scat- 
tered at an angle 1?, Ay^p^d) is the analyzing power of 
the p + CH2 reaction, and ip is the azimuthal scatter- 
ing angle. The additional terms ci , si , C2 , S2 , . . . represent 
false or instrumental asymmetries caused by non- uniform 
acceptance or efficiency, and possible (/s-dependent re- 
construction errors. These terms depend on p, and 
the incident proton trajectory, on which the geomet- 
ric acceptance depends. Normalized angular distribu- 
tions n± = N^{ip)/N^ can be defined for each helic- 
ity state. The helicity-sum distribution + n_ can- 
cels the helicity-dependent asymmetries corresponding 
to the transferred polarization, providing access to the 
false asymmetries, while the helicity-difference distribu- 
tion n_|_ — n_ cancels the helicity-independent false asym- 
metries, providing access to the physical asymmetries. 

False asymmetry effects are strongly suppressed in the 
extraction of the transferred polarization components by 
the rapid (30 Hz) beam helicity reversal, which cancels 
the false asymmetry contribution (to first order) and also 
cancels slow variations of luminosity and detection effi- 
ciency, resulting in the same effective integrated luminos- 
ity for each beam helicity state. Since the elastic scatter- 
ing cross section on an unpolarized proton target is in- 
dependent of electron helicity, equal numbers of protons 
incident on CH2 are detected for positive and negative 
beam helicities. In the GEp-II experiment, the numbers 



* In this context, "transverse" means orthogonal to the incident 

proton's momentum direction. 
^ In the assumed coordinate system, the 2 axis is along the incident 

proton momentum, while the x and y axes describe the transverse 

coordinates in relation to the proton trajectory and the detector 

coordinate system, as described in the text. 
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FIG. 6. (color online) Focal-plane helicity-diflterence asymme- 
try n+-n^ = {NMns/2) [N+{ip)/N+ - N-{ip)/No], where 
Nbins is the number of <p bins and N^{ip), N^^ are defined as 
in equation Q, for the three highest points from GEp-II. 
Curves are fits to the data. See text for details. 



of events in each helicity state were always found to be 
equal within statistical uncertainties at the 10"'' level. In 
a polarization transfer measurement, equal integrated lu- 
minosities for each beam helicity are not strictly required 
to robustly separate the physical from the instrumental 
asymmetries, since the angular distribution can be nor- 
malized to the number of incident protons for each he- 
licity state. Nonetheless, having equal numbers of events 
in each helicity state maximizes the statistical precision 
of the measured asymmetry while minimizing the sys- 
tematic uncertainty in its extraction. The false asymme- 
try coefficients determined from Fourier analysis of the 
helicity-sum distribution can be used to correct the resid- 
ual second-order effect of the false asymmetry, which is 
small compared to other uncertainties in the data of this 
experiment and therefore neglected (see section IIlIB4[) . 
Figure [S] shows the helicity-difference asymmetry ?i+ — 
for the three highest points from GEp-II, inte- 
grated over the range of polar angles with non-zero 
analyzing power. The data were fitted with n+ — n_ = 
a cos (/3 -I- 6 sin tp, with a resulting x^/ndi of 0.90, 0.53 and 
0.92 for Q2 = 4.0, 4.8 and 5.6 GeV^, respectively. At each 



Q^, the asymmetry exhibits a clear sinusoidal behavior, 
with a large cos ip amplitude proportional to P^^^ and 
a smaller sunp amplitude proportional to Py^^. There 
is no evidence in the data for a constant offset or the 
presence of higher harmonics, judging from the good 
of the fit with only cos ip and sin Lp termfl The ampli- 
tude of the asymmetry is proportional to the product 
of the weighted-average analyzing power and the magni- 
tude of the proton polarization, while the phase of the 
asymmetry is determined by the ratio Py^^/P^^^ of 
the proton's transverse polarization components at the 
focal plane. 



2. Spin Precession 

The asymmetry measured by the FPP is determined 
by the proton's transverse polarization after undergoing 
spin precession in the magnets of the HRS. To extract the 
transferred polarization components at the target corre- 
sponding to equations ([1]) requires accurate knowledge of 
the spin transport properties of the HRS. It is worth not- 
ing that without spin precession in magnetic spectrom- 
eters, a common feature of the GEp-I, GEp-II, GEp-III 
and GEp-27 experiments, proton polarimetry based on 
nuclear scattering would not work, since the spin-orbit 
coupling responsible for the azimuthal asymmetry is in- 
sensitive to the proton's longitudinal polarization, which 
can only be measured by rotating the longitudinal com- 
ponent into a transverse component. 

The precession of the spin of particles moving relativis- 
tically in a magnetic field is governed by the Thomas- 
B.M.T. equation [s^. The dominant precession effect 
in all of the aforementioned experiments is caused by the 
large vertical bend of the proton trajectory in the dipoles 
of the magnetic spectrometers. In first approximation, 
the proton spin precesses in the dispersive (vertical) plane 
by an angle x = li^p^bend relative to the proton trajec- 
tory, where 7^ = 1 -I- p^/M^ is the proton's relativistic 
boost factor, Kp is the proton's anomalous magnetic mo- 
ment, and Obend is the vertical trajectory bend angle. In 
this idealized approximation, the proton spin does not 
precess in the horizontal plane. The sensitivity of the 
FPP asymmetry to Pg is maximized when |sinx| = 1. 
The central values of x for the four kinematic settings of 
GEp-II are given in Table HI 

Because the central value of x is close to 360° at 
= 5.6 GeV^ and the range of x accepted by the HRSL 
is roughly 285° < x < 390°, the dominant cosp am- 
plitude of the focal plane asymmetry, which is roughly 
proportional to P^sinx, is reduced when averaged over 
the full X acceptance, as in the bottom panel of Figure 



Fits with Fourier modes up to and a constant term found 
that the coefficients of all terms other than cos ip and sin were 
zero within statistical uncertainties. 
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[6l However, the adverse impact of the unfavorable pre- 
cession angle on the precision of the data is mitigated 
by the large x acceptance of the HRS and the fact that 
Pi is quite large for the kinematics in question. The x~ 
dependence of the asymmetry is accounted for by the 
weighting of events in the unbinned maximum-likelihood 
analysis described below, which optimizes the statistical 
precision of the extraction without explicitly removing 
events near x = 360°. Moreover, the x and accep- 
tances of the HRSL are only weakly correlated, so that 
the range of contributing to the determination of R 
is not strongly affected. 

The presence of quadrupole magnets complicates the 
spin transport calculation by introducing precession in 
the horizontal (non-dispersive) plane, which mixes Pt 
and . The trajectory bend angle in the non-dispersive 
plane is zero for the spectrometer central ray, but non- 
zero for trajectories with angular and/or spatial devia- 
tions from the HRS optical axis. Because of the strong 
in-plane angle {9p) dependence of the cross section, the 
acceptance-averaged horizontal precession angle is gen- 
erally significantly non-zero. The quadrupole effects are 
qualitatively characterized by the non-dispersive preces- 
sion angle X4> = It^p't'bend, where 4>hend is the total tra- 
jectory bend angle in the non-dispersive plane. 

The spin transport calculation for the final analysis 
was performed using COSY (sJl, a differential algebra- 
based software library for charged particle optics and 
other applications. Since each proton trajectory through 
the HRS magnets is unique, the spin transport matrix 
must be calculated for each event. Rather than per- 
form a computationally expensive numerical integration 
of the BMT equation for each proton trajectory, a poly- 
nomial expansion of the forward spin transport matrix 
up to fifth order in the proton trajectory angles, vertex 
coordinates and momentum was fitted to a sample of 
random test trajectories that were propagated through a 
detailed layout of the HRS magnetic elements including 
fringe fields. The coefficients of this polynomial expan- 
sion were then used to calculate the spin rotation matrix 
for each event. Unlike the optics matrices used for parti- 
cle transport, which are independent of the HRS central 
momentum setting due to the fixed central bend angle, 
the spin transport matrix depends on the central momen- 
tum setting because the precession frequency relative to 
the proton trajectory is proportional to 7. Therefore, the 
fitting procedure for the COSY matrices had to be car- 
ried out separately for each . The Taylor expansion 
of the matrix elements in powers of the small deviations 
from the central ray within the acceptance of the HRSL 
converges quite rapidly to an accuracy better than the 
spectrometer resolution. 

Several coordinate rotations are involved in the cal- 
culation of the spin transport matrix elements for each 
event. First, the reaction plane coordinate system defines 
Pt and Pf. Pt is directed along the recoiling proton's mo- 
mentum and Pt is transverse to the proton momentum 
but parallel to the scattering plane, in the direction of de- 



creasing 9 p. A rotation is applied from the reaction plane 
to the fixed transport coordinate system in which the z- 
axis is along the HRS optical axis, the x-axis points along 
the dispersive plane in the direction of increasing parti- 
cle momentum (vertically downward), and the y-axis is 
chosen as y = i x 2; so that (a;, y, z) forms a right-handed 
Cartesian coordinate system. The COSY calculations are 
performed in this fixed coordinate system. After apply- 
ing the COSY rotation, which transports the spin from 
the target to the focal plane in transport coordinates, a 
final rotation is applied to express the rotated spin vec- 
tor in the comoving coordinates of the proton trajectory 
at the focal plane, in which the z axis is along the pro- 
ton momentum, the y axis is chosen perpendicular to the 
proton momentum and parallel to the yz plane of the 
transport coordinate system, and the x axis is chosen as 
X — y X z. The definition of the x and y axes of the co- 
moving coordinate system at the focal plane is arbitrary 
as long as it is applied consistently with the other coordi- 
nate systems involved. In the original analysis of GEp-II 
0, , the azimuthal FPP scattering angle (p was mea- 
sured counterclockwise from the y axis toward the x axis, 
as viewed along the z axis. This convention is also used 
in the present work, but it is worth noting that a different 
convention was used in the analysis of the GEp-III (2^ 
and GEp-27 [2^ experiments, in which was measured 
clockwise from the x axis toward the y axis. 

The observables Pt, Pi, and R were extracted from 
the data using an unbinned maximum-likelihood method. 
Up to an overall normalization constant independent of 
Pt and Pi, the likelihood function is given by 



C{Pt,Pi. 



n ^ i+Ao(^.)+ 

i=l L 

/^.Fe4'' {iS^}Pt + S^jPi) COS^. 



(s'^}Pt + s':^}Pi)s\np, 



(5) 



where Aq represents the sum of all false asymmetry terms, 
hi and P^ are the beam helicity and polarization, re- 
spectively. Ay is the analyzing power, and the Sjf! with 
j — x,y and k = t,i are the spin transport matrix el- 
ements. The values of Pt and Pi extracted by maxi- 
mizing the likelihood function ([S]) correspond to those of 
equations ([T]) in the case Pe = 1, i.e., the beam is 100% 
polarized. Converting the product over all events into 
a sum by taking the logarithm and keeping only terms 
up to second order in the Taylor-expansiorQ of the loga 



rithm (ln(l -|- x) 



^ /2 + Ox^, where x corresponds 



to the asymmetry) reduces the coupled, nonlinear system 



^ The maximum truncation error in the expansion of the loga- 
rithm for X = 0.1, an upper hmit corresponding to the largest 
i?-dependent asymmetries observed in the GEp-II data, is ap- 
proximately 0.3% (relative). 
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of partial differential equations to a linear system of alge- 
braic equations for the polarization transfer components: 



Pt 
Pi 



A«(l-A«) J'(') 



in which a sum over all events (X^i^i""*) implied, and 
the coefficients At and Xi are defined for the j*'' event as 

a|'' = h,P,A'^^ (4;' cos(^, - simp, 

Af = /;.Pe^,|,'^ [si') cos - sin . (7) 

Equation ^ can be written as a matrix equation 
MP = b, where M is the 2x2 matrix of sums multiply- 
ing the vector P of polarization transfer components, and 
b is the vector of sums on the right-hand-side of © • The 
solution of this equation is P = M~^h, and the standard 
statistical variances in Pt and Pi are obtained from the 
diagonal elements of the covariance matrix M~^. The 
corresponding statistical error in i? = ^pG^/G^j is ob- 
tained by appropriate error propagation through equa- 
tion ll]). The kinematic factor in equation Q is calcu- 
lated for each event from the reconstructed kinematics, 
and is averaged over all events in the calculation of R. 
Since the reconstruction of the kinematics is not unique 
and can be fixed by choosing any two of E^, E'^, 9^, Pp 
and 9p, the choice was made to use the quantities mea- 
sured with the highest precision, namely Pp and Eg, to 
calcul ate and e for each event. The kinematic fac- 
tor •\/r(l -|- e)/2e is known to a much better accuracy 
than the statistical and systematic accuracy of Pt / Pe and 
therefore makes a negligible contribution to the total un- 
certainty. 

It is worth remarking that "bin centering" effects due 
to the finite and e acceptance within each data point 
are essentially negligible, since the acceptance is 
small compared to the magnitude of Q^. The differ- 
ence between the average value of the kinematic factor 
yr(l + e)/2e and its value calculated at the average 
is negligible compared to the uncertainty in the ratio 
Pt/Pf- Furthermore, both the observed and expecte^l 
variations of Pt , Pe and R within the acceptance of each 
data point are small compared to their statistical un- 
certainties. Therefore, all data from each point are 
combined into a single result quoted at the average Q^. 

The forward spin transport matrix depends on all pa- 
rameters of the scattered proton trajectory before it en- 
ters the HRSL. Since the expected variation of R within 
the acceptance of each data point is small, any anoma- 
lous dependence of the extracted R on the reconstructed 
proton trajectory parameters is a signature of problems 
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FIG. 7. Dependence of extracted R = fipG^/G^j values on 
the reconstructed proton trajectory parameters at = 4.8 
GeV^: dtgt {4>tgt) is the proton trajectory angle relative to the 
HRS optical axis in the dispersive (non-dispersive) plane, ytgt 
is the position of the interaction vertex in spectrometer coor- 
dinates (see text for details), and 5 is the percentage deviation 
of the proton momentum from the central HRS momentum 
setting. In each panel, the data are integrated over the other 
three variables. 



with the spin transport calculation. Conversely, the ab- 
sence of anomalous dependence serves as a powerful data 
quality check. Figure [7] shows the dependence of R at 
= 4.8 GeV^, extracted using equation ([S]), on all four 
proton trajectory parameters that enter the spin trans- 
port calculation. These include the trajectory angles 
Stgt — ta.Yi~^ [dx / dz) and 4>tgt — ta.n~^{dy/dz) relative 
to the HRS optical axis, the vertex coordinate ytgt, de- 
fined as the horizontal position of the intersection of the 
proton trajectory with the plane normal to the HRS opti- 
cal axis containing the origir^, and S = 100 x {p — po)/PO: 
the percentage deviation of the measured proton momen- 
tum from the HRS central momentum setting. There is 
no evidence for a dependence of R on any of the variables 
involved in the precession calculation, indicating the ex- 
cellent quality of the COSY model. Linear and quadratic 
fits to the individual dependencies were also performed, 
and all non-constant terms included in the fits were found 
to be consistent with zero. 



Expected variations are based on the best current knowledge of 
the dependence of G'^/G^f. 



^ Assuming that the HRSL points at the origin of Hall A, ytgt is 
related to the position Zytx of the interaction point along the 
beamline by ytgt = —z^tx {sinQp + cosQptaniptgt), where 0p 
is the HRS central angle (given as Op in Table |T]|. 
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FIG. 8. Dependence of the form factor ratio R on the FPP 
polar scattering angle i3 for the three highest values of 
GEp-II. The constant behavior of R confirms the cancellation 
of Ay in the ratio Pt/Pi- 



3. Analyzing Power Calibration 

The p + CH2 analyzing power relating the size of the 
measured asymmetry to the proton polarization depends 
on the initial proton momentum and the scattering angle 
"d. Given the relatively small momentum acceptance of 
the HRS, the p-dependence of Ay within the acceptance 
of each point is much weaker than the very strong 1} 
dependence, and can be neglected as a first approxima- 
tion. Dedicated measurements of Ay (2^ at and above 
the momentum range of the GEp-11 experiment were per- 
formed prior to the GEp-111 experiment. However, pre- 
cise independent knowledge of Ay is not required in the 
analysis because of the self-calibrating nature of elastic 
ep scattering, explained below. 

Provided that the effective d acceptance is (p- 
independent, the analyzing power cancels in the ratio 
Pt/Pt from which the form factor ratio R is extracted, 
implying that the result for R is independent of Ay . Uni- 
form 1? acceptance is guaranteed by applying a "cone 
test" in the selection of FPP events, which requires that 
the projection to the rearmost FPP detector plane of a 
track originating at the reconstructed P-I-CH2 scattering 
vertex Zdose at a polar angle "i? falls within the active 
detector area for all azimuthal angles (p. Moreover, the 
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FIG. 9. (color online) Top panel: extracted analyzing power 
as a function of p sin where p is the proton momentum inci- 
dent on CH2 (corrected event-by-event for energy loss in CH2 
up to the reconstructed scattering vertex), using the Pe values 
of Table [II for all four values of the GEp-II experiment. 
Gurves are fits to the data (see text for details). Bottom 
panel: maximum analyzing power vs. 1/po in GeV~^, where 
po is the central proton momentum, for the four points. 
Error bars in Ay and A™"^ values are statistical only. See 
Supplemental Material at [S^ for data tables with numerical 
Ay and A™""" results. 



cancellation can be verified by binning the results in 
and checking the constancy of i? ex Pt/P^ as a function 
of d. Figure [5] shows the dependence of R for the three 
highest points of GEp-11. At each Q^, a constant fit 
to the data gives a good and no systematic trends are 
observed. 

The fact that Pt and Pe depend only on R and kine- 
matic factors implies that the product PeAy can be ex- 
tracted by comparing the measured asymmetries PgAyPt 
and PeAyPe to the values of Pt and Pe obtained from 
equation ([1]). Combined with the measurements of P^ to 
within an overall accuracy of ±3% by MoUer and Comp- 
ton polarimetery, Ay was directly extracted from the data 
of this experiment. The p and dependences of Ay thus 
obtained were then used in equation ^ to improve the 
statistical precision of the form factor ratio extraction by 
weighting events according to their analyzing power. 
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Vt (GeV) 


a 






3.5 


0.030 ± 0.008 


0.89 ± 0.06 


1.19 ±0.10 




4.0 


0.031 ± 0.003 


0.88 ± 0.03 


0.95 ± 0.04 




4.8 


0.029 ± 0.005 


1.02 ±0.04 


1.03 ±0.05 




5.6 


0.038 ±0.011 


1.14 ±0.09 


1.12 ±0.11 


Q 


(GeV ) 


7, 



X /n.d.t. 


A max 




3.5 


3.51 ±0.15 


1.23 


0.149 ±0.005 




4.0 


3.28 ±0.06 


1.22 


0.123 ±0.003 




4.8 


3.44 ± 0.06 


1.30 


0.109 ±0.002 




5.6 


3.67 ±0.15 


2.02 


0.099 ± 0.004 



TABLE II. Ay fit results. Parametrization is Ayipr) = 

(pT -pT)"e-''(^^-^?■)^ The uncertainty in A™" was cal- 
culated from the full covariance matrix of the fit result. 



Figure [S] shows the measured Ay as a function of the 
"transverse momentum" px = psint? for each point, 
where p is the incident proton momentum corrected for 
energy loss in CH2 up to the reconstructed scattering ver- 
tex, illustrating the approximate scaling of the angular 
distribution of Ay with momentum. The results shown 
in Figure [S] are in fairly good agreement with the unpub- 
lished results from the original analysis in [295 ^ despite us- 
ing the more restrictive elastic event selection cuts of the 
present work. This is due in part to the fact that the sen- 
sitivity of Pf, from which Ay is primarily determined, to 
r = G^^/G\j is rather weak (see equation ([1])). Nonethe- 
less, for the three highest points, the improved sup- 
pression of the background in this analysis leads to a 
slight systematic increase in Ay, since the asymmetry of 
the background included in the original analysis partially 
cancels that of the signal. Ay rises rapidly from zero in 
the region dominated by Coulomb scattering to a maxi- 
mum at pt ~ 0.3 GeV and then tapers off to nearly zero 
beyond about 1.5 GeV. The measured angular distribu- 
tion at each was fitted using a simple parametrization 
Ay{pT) = (pT - p^)"'e-''^P^-PT'>\ where p^, a, b and 
/3 are adjustable parameters. This parametrization in- 
corporates the main features of the angular distribution 
with sensible limiting behavior and is sufficiently flexible 
to give a good description of the data. The fit results 
for each are given in Table [TTl The quality of the 
fit was improved by including the zero offset the 
data seem to prefer a vanishing Ay at finite p^ « 0.03 
GeV, independent of Q^. For pT < Pj', Ay — was 
assumed. The results for the exponents a and /3 are 
essentially compatible with the product of a linear rise 
and an exponential decay. An alternate parametrization 
which fixes a = 1 and /? = 1 and adds an overall nor- 
malization constant as a free parameter in addition to 
the slope parameter b does not describe the data as well 
as the chosen parametrization in which a and (3 are free 
parameters but the overall normalization is fixed. The 
amplitude of the measured Ay distribution, as measured 
by its maximum value, scales approximately with 1/p, 
as shown in the bottom panel of Figure IH) Notably, the 
intercept of the linear fit to the 1/p dependence of yl™"^ 



is compatible with zero, suggesting that the analyzing 
power for p + CH2 scattering vanishes for asymptotically 
large proton momenta, rather than crossing zero at a fi- 
nite momentum. The fitted curves shown in Figure [3] 
were used to describe Ay{pT) in the analysis. 

The observed proportionality of Ay to 1/p allows the 
momentum dependence of Ay to be accounted for in the 
analysis by simply scaling its value for each event by a 
factor po/p, where po is the central proton momentum 
and p is the proton momentum for the event in ques- 
tioifl. This is because the fitted Ay{pT) curve, which 
is averaged over the ±5% momentum bite of the HRS 
at each Q^, essentially gives Ay{po,PT), where po is the 
central momentum. Assuming that the 1/p slope of Ay 
is the same at any px] i.e., assuming a factorized form 
^v{PtPt) = C{pt)/p, the ratio of Ay{p,pT) to its known 
value Ay{po,pT) at a reference momentum po is given by 
Po/p, regardless of C{pt)- While the observed shape of 
the Pt dependence of Ay is approximately momentum- 
independent for the three higher points, the pr de- 
pendence of Ay at = 3.5 GeV^ is slightly different, 
with a larger maximum value than suggested by a lin- 
ear extrapolation from the higher-Q^ data and a faster 
falloff at large px- A plausible, but unproven explana- 
tion for the difference in behavior is that the thicker 100 
cm analyzer used for the three highest-Q^ measurements 
smears out the pt distribution of both the efficiency and 
the analyzing power of the FPP relative to the thinner 
58 cm analyzer used for the measurement at = 3.5 
GeV^. This observation does not, however, invalidate 
the po/p scaling of Ay in the analysis, because the data 
from the three higher- points, as well as data from 
other experiments [J, [2^, show that the 1/p scaling is 
respected for any given FPP configuration, though the 
details of Ay{pT) may differ slightly between different 
configurations. In any case, the value of Ay assigned in 
the analysis is never changed by more than ±5% for any 
individual event, so the actual effect of this prescription 
on the relative weighting of events is rather small. 

The description of Ay {p, 1?) in the present reanalysis 
differs slightly from that of the original analysis. In this 
reanalysis, Ay{p, is assigned to each event based on the 
smooth parametrization of Ay{pT) shown in the curves 
of Figure [51 which describe the data very well, and an 
overall 1/p scaling. The original analysis, on the other 
hand, neglected the momentum dependence of Ay and 
assigned Ay{'d) to each event based on the calibration 
results in discrete "i? bins. Since Ay cancels in the ra- 
tio Pt/Pg, its description only matters to the extent that 
it optimizes the statistical precision of the extraction. 
Different descriptions of Ay{p,'d) correspond to different 
event weights in the analysis, leading to slight differences 



For this purpose, the central momentum po was corrected for 
energy loss in half the thickness of CH2 , while the momentum p 
for the event in question was corrected for energy loss up to the 
reconstructed scattering vertex. 
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in the results for Pt, Pi and R reflecting statistical fluc- 
tuations of the data as a function of p and d. While 
these differences are always well within the statistical 
uncertainty of the combined data, better descriptions of 
Ay{p,d) naturally lead to better overall results. 

4. False Asymmetries 

Consistent with the original analysis, no false asym- 
metry corrections were applied in the present work; i.e., 
Ao = was assumed in equations ([5]) and ^ . "Weighted 
sum" estimators, as defined in [s^, can be constructed for 
the focal plane asymmetries Ay^^ = —PeAyPy^^ and 
A^^^ = PeAyPj^^^ , equivalent to equation (jH) in the 
absence of precession effects. Including false asymmetry 
terms up to 2ip, it can be shown that the weighted-sum 
estimators A^^^ and Ay^^ for the focal plane asymme- 
tries are given to second order in the false and physical 
asymmetry terms by 

ir^ = Ar^(i-|)-<^^| 

A^^^^-'-^Ar^ + A^^^{l + (8) 

where C2 and S2 are the false asymmetries as in equa- 
tion (H]). Only the 2(p Fourier moments of the false 
asymmetry contribute at this order. The cos(2(^) false 
asymmetry moment induces a "diagonal" correction to 
each physical asymmetry term proportional to the asym- 
metry itself, while the sin(2<p) false asymmetry moment 
induces an "off-diagonal" correction to A^^^ (Ay^^) 
proportional to A^^^ {A^^^). 

Fourier analysis of the helicity sum distribution n+ + 
n_ showed that the acceptance-averaged magnitude of 
C2 and S2 did not exceed 2.5 x 10~^ at any Q^, and nei- 
ther term exceeded 1% at any d within the useful range. 
The possible effect of C2 on the "diagonal" terms is there- 
fore at the 10~^ (relative) level, while the "off-diagonal" 
correction is at the 10~^ level (absolute) for the small 
Ay^^ term, and even smaller for the larger A^^^ term. 
Compared to both the size and statistical uncertainty in 
the asymmetries (see Fig. [6]), and the systematic uncer- 
tainties in Pt and Pi resulting from the spin transport 
calculation, such corrections are completely negligible. 
This is in contrast to the GEp-III and GEp-27 analyses, 
in which a sizeable cos(2(p) false asymmetry in the Hall 
C FPP induced a correction that, while small, made a 
non-negligible contribution to the total systematic un- 
certainty. 

C. Background Estimation and Subtraction 

From Figure [5l two qualitative features of the data 
are obvious. First, the non-elastic background before ap- 
plying two-body correlation cuts is substantial. Second, 
examination of the Ax and Ay spectra before and after 



applying the 5p cut reveals that the 5p cut provides signif- 
icant additional background suppression power relative 
to Aa; and Ay cuts alone, with minimal reduction of the 
elastic peak strength, implying that events outside the 
Sp cut are background-dominated, even after calorimeter 
cuts. 

As alluded to in sections [II Bl and IIII Al the non-elastic 
background for the measurements using a calorimeter for 
electron detection consists predominantly of two reac- 
tions, quasi-elastic Al(e, e'p) scattering in the cryocell en- 
trance and exit windows, and tt^ production initiated by 
the flux of real Bremsstrahlung photons radiated along 
the target material (photoproduction) as well as virtual 
photons present in the electron beam independent of tar- 
get thickness (electroproduction). Due to the kinematic 
acceptance of the experiment and the dependence of 
the respective cross sections, the contribution oin'^p elec- 
troproduction is mostly limited to "quasi-real" photons; 
i.e., ss 0, and is practically indistinguishable from 
real photoproduction. By detecting both scattered par- 
ticles in coincidence, the two-body ep — >■ ep kinemat- 
ics are overdetermined, providing for a clean selection of 
elastic events and a direct determination of the remain- 
ing background from the data, with no external inputs, 
using the sideband-fitting method described in section 
IIII C 21 below. The main disadvantage of this approach 
to background estimation is that it makes no reference 
to the underlying physics of the signal and background. 
For this reason, a Monte Carlo simulation of the experi- 
ment was carried out to confirm the conclusions regarding 
backgrounds obtained directly from the data. However, 
the results of the simulation were not used in any way as 
input to the final analysis. 



1. Monte Carlo Simulation 

The simulation code is the same as that used in the 
data analysis of which already includes a realistic 
model of the HRSL. Modifications of the code used in 
the analysis of Ref. [Ej to reproduce non-Gaussian tails 
of the HRS resolution, caused by multiple scattering and 
other effects, were not included here. The only signifi- 
cant addition to the code was a description of the accep- 
tance and resolution of the GEp-II calorimeter. Because 
the 15 X 15 cm^ cell size of the GEp-II calorimeter is 
large compared to the Moliere radius of lead-glass, coor- 
dinate reconstruction essentially consists of assigning the 
shower coordinates to the center of the cell with maxi- 
mum energy deposition. Furthermore, the discriminator 
threshold applied to form the timing signal was roughly 
20% of the elastically scattered electron energy, meaning 
that signals below this amplitude would be rejected in 
software by the timing cut. The electron energy and co- 
ordinates were thus defined by the signal in a single block 
in the overwhelming majority (> 90%) of elastic events. 
Physics ingredients of the simulation include cross sec- 
tion models for ^H(e, e'p), Al(e, e'p) and ^H(7, n^p) reac- 
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FIG. 10. (color online) Contributions to the Sp distribution at 
= 4.8 GeV^ estimated from the Monte Carlo simulation. 
Monte Carlo distributions are shown for elastic ep (empty 
circles), jp — i- n^'p (empty triangles), and quasi-elastic (e, e'p) 
in the target windows (empty inverted triangles). The sum of 
all Monte Carlo contributions (empty squares) is compared to 
the data of Figure[S{c) (solid inverted triangles). Monte Carlo 
and data distributions are obtained after applying Aa; and Ay 
(calorimeter) cuts. Black vertical lines show the Sp cut region 
of the final analysis. Uncertainties shown are statistical only. 
See text for details. 



tions, a realistic calculation of the Bremsstrahlung flux 
for 7r° photoproduction, and event-by-event radiative cor- 
rections to the (e, e'p) cross sections following the ap- 
proach of [37], providing for a rigorous deconvolution of 
the signal and background contributions to the Ax, Ay, 
and 6p distributions for arbitrary cuts. Another reaction 
that can contribute to the background is Real Compton 
Scattering — > jp (RCS), whose end-point kinemat- 
ics are identical to ep — t- ep. However, the cross section 
for this reaction is g enerally much smaller than for tt^ 
photoproduction [38l. [39| , and was neglected. 

Figure [TO] shows the simulated dp distribution in the 
vicinity of the elastic peak for each reaction considered, 
after applying Aa; and Ay cuts. As described below, the 
simulated target window yield was normalized to match 
the window yield obtained from the data in the super- 
elastic {Sp < 0) region. Then, the overall normaliza- 
tion constants for ir'^p and elastic ep events were fitted 
simultaneously to minimize the statistics-weighted sum 
of squared differences between the data and the sum of 
Monte Carlo yields. The agreement between data and 
Monte Carlo is good, but not perfect, primarily because 
non-Gaussian tails are not included in the simulated Sp 
resolution. Nonetheless, the Sp distribution after cuts is 
described to within ^ 20% in the relevant Sp range, with 



the exception of disagreements of up to ~ 40% in the Sp 
region from 20-40 MeV just above the elastic peak, which 
is rather sensitive to non-Gaussian tails and the details 
of the Bremsstrahlung spectrum and the tt^ production 
cross section near end-point. Since the purpose of the 
simulation was to provide a qualitative illustration of the 
physics of the signal and the background, and since the 
background contamination and its polarization were de- 
termined directly from the data for the final analysis, no 
additional fine-tuning of the simulation was attempted. 

Two key features of the simulation results deserve spe- 
cial emphasis. First, the contribution of the ep radiative 
tail in the inelastic region falls off too quickly to describe 
the observed tail of the data. This is a consequence of 
the Aa; cut, with Aa; calculated using the GEp-II method 
2^3. The background fraction exceeds 80% above 50 



MeV and 90% above 75 MeV. The ep yield falls below 
the tt'^p yield at ~40 MeV and becomes negligible above 
~ 120 MeV, confirming the conclusion that the inelas- 
tic region of the Sp distribution is dominated by the ir'^p 
background rather than the ep radiative tail. Second, 
the target window contribution is vanishingly small com- 
pared to the elastic and ir^p contributions in the entire Sp 
range of interest. More specifically, in the region below tt'^ 
threshold, the window contribution is the dominant com- 
ponent of the background, but is too small relative to the 
elastic yield to affect the measured asymmetry, while in 
the region where the contamination is sufficiently large 
to affect the asymmetry, the ir^ contribution is domi- 
nant. Moreover, the proton recoil polarization in quasi- 
elastic Al(e, e'p) scattering at high should be similar, 
in principle, to that in elastic ep — > ep, since the for- 
mer process is simply the latter process embedded in a 
nucleus, whereas the spin structure of 7p — > Tr'^p can be 
(and is) dramatically different. 

The only kinematically allowed reactions producing 
protons in the super-elastic region are quasi-elastic 
Al(e, e'p) and other reactions occurring on the Al nu- 
clei in the cryocell windows, in which the initial Fermi 
motion of the struck proton can lead to proton knockout 
with pp > Pp{9p). However, a significant fraction of the 
yield in the super-elastic region actually comes from hy- 
drogen, because the combined thickness of the entrance 
and exit windows of the Hall A cryotarget [l^] in g cm~^ 
is only about 4% of the liquid hydrogen thickness, and the 
non-Gaussian tails of the Sp resolution smear a fraction 
of hydrogen events into the unphysical Sp region. The 
reconstructed vertex distribution in this region exhibits 
narrow peaks at the window locations and a smooth hy- 
drogen background extending over the full target length. 
To estimate the yield from the target windows, the vertex 
z distribution was plotted as a function of Sp in the super- 
elastic region for events failing the Aa; and Ay cuts, in 
order to enhance the very small window "signal" relative 



The Ax cut suppresses the ep radiative tail even more strongly 
when Ax is calculated using the GEp-III method. 
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to the large hydrogen elastic "background". For each 
of six 5p bins in —180 < 6p (MeV) < 0, a polynomial 
fit to the smooth hydrogen background was subtracted 
from the vertex z distribution, leaving only the window 
peaks. For each window, the simulated Sp distribution 
with identical cuts applied was normalized to match the 
background-subtracted window yield obtained from the 
data. The resulting normalization factor was then ap- 
plied to the simulated Sp distribution of window events 
passing the Ax and Ay cuts, leading to the contribution 
shown in Fig. (TU] 

Given the vertex resolution of the HRS, a vertex cut 
chosen to exclude the windows at the 3a level can further 
suppress the very small window background, at the ex- 
pense of a ~20% reduction in elastic ep statistics. How- 
ever, the aforementioned analysis of the window yield 
suggests that even when the full target length is included, 
the fraction of the total yield from the windows is negli- 
gible after all cuts are applied, making additional ver- 
tex cuts unnecessary. This conclusion is further sup- 
ported by comparing the 6p distributions with and with- 
out such a vertex cut, and by comparing the Sp spectra 
for the > 4.0 GeV^ settings to the Sp spectrum of the 
= 3.5 GeV^ setting, for which the precise measure- 
ment of the electron kinematics with a magnetic spec- 
trometer provides an essentially background-free selec- 
tion of elastic events, as discussed in Section lll A2I Based 
on these considerations, the window contamination was 
deemed negligible, and the study of the background con- 
tamination focused mainly on the inelastic (Sp > 0) re- 
gion. 

The background subtraction procedure used for the fi- 
nal analysis is agnostic regarding the reaction mechanism 
responsible for the contamination, with the caveat that 
the conclusion of negligible window contamination is used 
to justify the assumption of constant background polar- 
ization, which reduces the statistical uncertainty in the 
background correction. In summary, the simulation pro- 
vides a qualitative description of the data that supports 
the conclusions of this analysis regarding backgrounds. 
Averaged over the final Sp cut region, the fractional 
background contamination obtained from the simulation 
agrees with that obtained directly from the data at a level 
similar to its systematic uncertainty, which is determined 
by the data. 



2. Sideband subtraction 

For the final analysis, the fractional background con- 
tamination in the sample of elastic ep events selected by 
a given set of cuts was estimated by fitting the tails of 
the Ax and Ay distributions on either side of the elastic 
peak and extrapolating into the peak region, as shown in 
Figures[5ja) and[5ljb). This approach to background esti- 
mation implies two assumptions. First, the contribution 
of elastic scattering to the tails of the Ax and Ay distri- 
butions is assumed to be negligible for values of Ax and 



Ay sufficiently far away from the elastic peak. Second, 
the background is assumed to have a smooth distribu- 
tion under the elastic peak, so that joining the tails with 
a smooth interpolating function is a good approximation 
to the true background shape. The first assumption can 
in principle be violated by the ep radiative tail and by 
non-Gaussian smearing effects in the HRS angle and mo- 
mentum reconstruction. Radiation redistributes elastic 
ep events away from the elastic peak toward negative Ax 
values, but does not markedly affect the Ay distribution 
of elastic events, since Ay reflects the extent to which the 
two detected particles are non-coplanar, and the copla- 
narity of outgoing particles is not strongly affected by 
radiation. Furthermore, the Sp cut suppresses the radia- 
tive tail of the Ax distribution. Non-Gaussian smearing 
effects do not contribute a significant fraction of events in 
the tails except when the background contribution is very 
small. The second assumption (smooth background dis- 
tribution) was confirmed by inspecting the correlations 
between Ax and Ay; i.e., by plotting Ax (Ay) for Ay 
(Ax) well outside the elastic peak. This assumption was 
also supported by the simulations described in section 
IIII C II Although the simulation does not include the 
contribution of random coincidences, the contamination 
of the data by random coincidences is negligible after 
timing and kinematic cuts. 

In the following discussion, the fractional background 
contamination / is defined as / = B/{S + B), where B 
is the number of background events and S is the num- 
ber of signal events; i.e., / is the ratio of the background 
yield to the total yield. The value of / and its system- 
atic uncertainty A/ were estimated using a conserva- 
tive approach involving a total of twelve different fits. 
The tails of the Ax and Ay distributions, obtained after 
applying all other cuts, were each fitted with Gaussian 
and polynomial background shapes, for three different 
sizes of the elastic peak region excluded from the fit (two 
spectra x two parametrizations x three sideband ranges 
= twelve fits). The average fit result was taken as the 
value of /, while the rms deviation of the fit result from 
the mean was taken as the systematic uncertainty A/. 
The variations among the different fit results reflect the 
level of agreement (or disagreement) among the differ- 
ent spectra, assumed background lineshapes, and regions 
excluded from the fit. 

A central conclusion of the present reanalysis is that 
the background was underestimated in the original anal- 
ysis. Using the polynomial sideband fitting method, the 
estimated average values of / for the cuts of the original 
analysis, in which no Sp cut was applied, are 1.6%, 2.8% 
and 5.3% for = 4.0, 4.8 and 5.6 GeV^, respectively. 
Compared to the estimates reported in j2^ for the origi- 
nal analysis, these estimates are higher by factors of 2.3, 
7.0 and 3.8, respectively. Even at the few percent level, 
neglected or underestimated inelastic contamination can 
have a non-negligible effect on the measured asymme- 
tries if the polarization of the background differs strongly 
enough from that of the signal, as in this case. 
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With the addition of the Sp cut, the present analy- 
sis maximally exploits the two-body kinematic correla- 
tions of both detected particles. In the inelastic region, 
tt'^ production dominates. In terms of Sp, the tt^ pro- 
duction "threshold" is very close to the elastic peak. 
When reconstructed assuming elastic scattering, protons 
from 7p — TT^p at the Bremsstrahlung end-point have 
5p = 7.4, 8.1 and 8.8 MeV for = 4.0, 4.8 and 5.6 
GeV^, respectively. When compared to the dp resolution 
of ~ 5 MeV, there is clearly substantial overlap of the 
TT^p kinematic phase space with the elastic peak, as in 
the example of Figure [TUl As increases at a given 
beam energy, the n'^p cross section becomes large com- 
pared to the ep cross section. 

The effect of underestimating the tt" background on 
the form factor ratio extraction is illustrated in FigurelTTl 
which shows Pt, Pe and / as a function of Sp, for events 
identified as elastic in the original analysis, at — 4.8 
GeV^. The data were divided into eight Sp bins, includ- 
ing six equal-statistics bins inside the cut region of Fig- 
urelSJc), where / is very small (—7.3 < Sp < 22.7 MeV), 
a seventh bin with a significant fraction of both signal 
and background (22.7 < Sp < 60 MeV), and an eighth 
bin dominated by background {Sp > 60 MeV). Because 
the Aa; and Ay distributions in the last Sp bin showed 
no obvious signature of an elastic peak, / = 1 was as- 
sumed for this bin, consistent with the simulation results 
shown in Figure [TUl Meaningful background estimation 
and subtraction were not possible for this bin. As Sp 
increases, the raw transferred polarization components 
P°'"* and P^°''" evolve from their roughly constant values 
in the signal-dominated region to values that are consis- 
tent with the background polarization components P™'^^ 
and Pp^K The (5p-integrated results for the background 
polarization, extracted from events rejected by the cuts 
of Figure [SI are plotted at an arbitrary Sp = 115 MeV 
for comparison. 

The background polarization components were ob- 
tained by applying anti-cuts twice as wide as the final 
elastic event selection cuts; i.e., Ax{Ay) was required to 
be at least 24(32) cm away from the midpoint between 
half-maxima of the peak. Events selected by this anti- 
cut are background-dominated and have negligible elas- 
tic contamination. In order to study the Sp dependence 
of no cut was applied to Sp in the extraction of 

the background polarization. No statistically significant 
Sp dependence of the background polarization was ob- 
served, consistent with dominance of the background by 
7r°p events. Therefore, P?""^' was assumed constant in 
the background subtraction procedure. 

In Figure [TTJ the signal polarization P/'(i — t,£) was 
obtained from P°'"* in the first seven bins using the sub- 
traction 
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FIG. 11. (color online) Top panel: 5p dependence of Pt and 
Pi for events selected by the Aa; and Ay cuts of Figure [S] 
Data are binned in 5p as described in the text and plotted 
at the average 5p value in each bin. Raw polarizations P°''" 
(filled squares) and P^''" (filled circles) approach the back- 
ground polarizations Pt""'' (filled triangle) and p™^' (empty 
triangle) at large Sp. Corrected values Pt' (empty squares) 
and Pg' (empty circles) are ofi^set in 5p for clarity. Dashed and 
solid horizontal lines are weighted averages of the corrected 
and raw data, respectively. Bottom panel: Sp dependence of 
the fractional background contamination /. Uncertainties in 
/ are systematics-dominated, while the uncertainties in the 
polarization components are statistics-dominated. See text 
for details. 



By comparing the weighted average of all uncorrected 
data in Figure [TT] to the weighted average of the six 



corrected data points inside the cut region, it is found 
that the background contamination of the sample with 
no Sp cut induces relative systematic shifts of |APf/Pf | = 
15.8% and \ APi/Pi\ = 2.4%. From Figure [IH it is clear 
that the tails of the Sp distribution outside the cut region 
of Figure [5l[c) contribute very little to the statistical pre- 
cision of the measurement of Pt /Pe while causing a large 
systematic effect. For the final analysis, rather than cor- 
recting the results bin- by-bin in Sp using equation ([SJ, 
as in Figure [TTJ the background fraction / and polariza- 
tion p?"'^' were included at the individual event level in 



18 



equation ([6]) by making the following replacements: 



1-A, 



A 



inel 



(10) 



where fi is the background contamination as a function of 
(5p(*^ and A -^^j , representing the background asymmetry, 
is given by 



A« , ^ cos^, - sin^,) pr^'- 



( 



S'^'^'' cos ipi - Sl"^ sin ) Pg 



(11) 



This method is functionally equivalent to correcting the 
results "after the fact" using equation It also simpli- 
fies the evaluation of systematic uncertainties associated 
with the background correction, which were obtained by 
varying /, p™«' and P^*"*^' within their uncertainties and 
observing the shift in R. 



D. Systematic Uncertainties 

As a result of the cancellation of the beam polarization 
and analyzing power in the ratio Pt/Pi and the cancel- 
lation of the FPP instrumental asymmetry by the beam 
helicity reversal, there are few significant sources of sys- 
tematic uncertainty in the results of this experiment (as 
is also the case in the GEp-I, GEp-III and GEp-27 exper- 
iments). The dominant source of systematic uncertainty 
is the spin transport calculation. Since the procedure 
for the evaluation of systematic uncertainties associated 
with this calculation is documented at length in Refs. 
[U, m, 113, im : only a brief summary of the studies and 
the conclusions is given here. 

The range of non-dispersive plane trajectory bend 
angles 4>bend accepted by the HRS is roughly ±60 
mrad, independent of momentum. The maximum ac- 
cepted range of the non-dispersive plane precession an- 
gle X0 = jKp4>bend is roughly ±30° at the highest 
of 5.6 GeV^. To first order in x<p^ the ratio Pt/Pe is 
given in terms of the focal plane ratio Py^^/P^^^ 

by PtjPg ~ X<i> ~ ^^"^xPy^^ I Px^^ ■ Because the non- 
dispersive plane precession mixes Pt and P^, the ratio 
is highly sensitive to uncertainties in (pbend- To first 
order, an uncertainty A(j)hend leads to an uncertainty 

AP ss ^^p-y/r(l + e)/2ej "fKp/S.(f)bend in the extracted 
form factor ratio. The error magnification factor multi- 
plying A(j)i,end grows as large as 33 at = 5.6 GeV^. To 
manage the systematic uncertainty due to the precession 
calculation, (fibend must be known to very high accuracy. 
On the other hand, since Obend only enters Pt / Pi through 
the factor of sinx multiplying Py^^/P^^^, and since 
the reconstruction of Obend involves relatively small de- 
viations about the 45° central bend angle, the accuracy 
of Pt /Pi is far less sensitive to systematic errors in Obend 
and Pp . 



The major sources of uncertainty in 4>bend are horizon- 
tal misalignments and rotations of the three quadrupoles 
relative to the HRSL optical axis defined by the dipole 
magnet. In order to control the uncertainty in (pbend to 
the highest possible accuracy, dedicated studies of the 
optical properties of HRSL in the non-dispersive plane 
were performed. Electrons were scattered from a thin 
carbon foil aligned with the HRSL optical axis, and a 
special "sieve-slit" collimator was installed in front of 
the entrance to HRSL before the first quadrupole mag- 
net. The sieve-slit collimator, part of the standard equip- 
ment of the HRSs, consists of a 5 mm thick stainless steel 
sheet with a pattern of 49 holes (7 x 7), spaced 25 mm 
apart vertically and 12.5 mm apart horizontally, used for 
optics calibrations '2?]. In the studies described here, 
electrons passing through the central sieve hole aligned 
with the HRS optical axis were selected. For a series of 
deliberate mistunings of the HRS quadrupoles relative 
to the nominal tune, the displacements in both position 
and angle of the image of the central sieve hole at the 
focal plane were observed. Combined with the known 
first-order HRS optics coefficients describing the effects 
of quadrupole misalignments and rotations, the informa- 
tion gained from these studies placed a much more strin- 
gent constraint on the misalignments than the nominal 
accuracy of the quadrupole positions. By reducing the 
uncertainty A(j)bend to ±0.3 mrad, the optical studies re- 
duced the systematic uncertainty in RaX — 5.6 GeV^, 
where the result is most sensitive to (pbend, to a level com- 
parable with other contributions. 

Additional model uncertainties in the precession calcu- 
lation due to the field layout in COSY are more difficult 
to quantify, but are typically smaller than the errors asso- 
ciated with the accuracy of the inputs to the calculation; 
i.e., the reconstructed proton kinematics. The COSY 
model uncertainties were estimated by performing the 
calculation in several different ways. For the final analy- 
sis, the proton trajectory angles, momentum and vertex 
coordinates, calculated using the standard HRS optics 
matrix tuned to calibration data as described in [27] , were 
used to calculate the forward spin transport matrix, as 
described in section UlI B 21 To estimate systematic un- 
certainties, the calculation was also performed using the 
same forward spin transport matrix, but the kinematics 
were reconstructed using an alternate set of optics ma- 
trix elements calculated by COSY. Finally, COSY was 
used to calculate the expansion of the reverse spin trans- 
port matrix, which was then inverted to obtain the for- 
ward matrix elements that enter the likelihood function 
of equation ([5]). A model systematic uncertainty was as- 
signed based on the variations in the results among the 
different methods, as described in [29j . 

Apart from the uncertainties associated with the non- 
elastic background, which were underestimated by the 
original analysis, the main additional source of uncer- 
tainty is the accuracy of the scattering angle reconstruc- 
tion in the FPP. Uncertainties associated with FPP re- 
construction were minimized by a software alignment 
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procedure using "straight-through" data obtained with 
the CH2 analyzers removed. The systematic errors in R 
due to the absolute accuracy in the determination of the 
beam energy {5E/E ^ 2 x 10^'') and the proton mo- 
mentum {Sp/p ^ 4 X 10^'*) [27], which mainly enter the 
ratio R through the kinematic factor /ip-y/r(l + e)/2e of 
equation ([T|), are negligible compared to the precession- 
related uncertainties. 

The updated systematic uncertainties associated with 
the background estimation and subtraction procedure are 
very small as a result of the added Sp cut, and are gener- 
ally at the 10~'^ level. The "Bckgr." uncertainty in Table 
mil was obtained by varying /, p™'^' and P™"^' within 
their uncertainties, which are systematics-dominated for 
/ and statistics-dominated for p™'='^ and observing the 
shift in R. The contributions from / and P™"^' are 
comparable, while the contribution from p™"^' is much 
smaller. 

The present analysis also examined the sensitivity of R 
to variations in elastic event selection cuts. The analysis 
was performed for various Ax, Ay and 6p cut widths, 
using both the GEp-II and GEp-III definitions of Ax and 
Ay (see section UlI Ap . The analysis was also performed 
using the original polygon cut, supplemented by the new 
Sp cut. For consistency of background corrections, the 
contamination was estimated separately for each case. 
The rms variation of R due to cut variations is quoted 
as the "Cuts" uncertainty of Table IIIII It is generally 
larger than the "Bckgr." uncertainty calculated using 
the final cuts, and reflects fluctuations among slightly 
different selections of events, not necessarily related to 
the background. It is, however, much smaller than the 
statistical uncertainty at each Q^. 

The present reanalysis of the GEp-II data is identi- 
cal to the original analysis in event reconstruction, spin 
transport calculations, and all cuts other than Aa;, Ay 
and Sp used to select elastic events. The only other mean- 
ingful difference between the present reanalysis and the 
original analysis is the improved description of the ana- 
lyzing power discussed in section IIII B 3[ which only af- 
fects the results through slight modification of the p and 
i?-dependent weighting of events. Therefore, aspects of 
systematic uncertainty analysis other than elastic event 
selection and background subtraction were not revisited. 
These aspects of the analysis are documented at length 
in [2|. 

Table IIIII shows all known contributions to the sys- 
tematic uncertainty in R at each Q^, including the po- 
lar {'&FPp) and azimuthal {(pppp) angle reconstruction 
in the FPP, the dispersive (Obend) and non-dispersive 
{4>bend) trajcctory bend angles, the COSY model uncer- 
tainty (COSY), the non-elastic background contribution 
(Bckgr.) and the cut sensitivity (Cuts). AU contribu- 
tions are added in quadrature to obtain the total system- 
atic uncertainty. Uncertainties due to FPP instrumental 
asymmetries are negligible as discussed in section IIIIBI 
In the final analysis of the GEp-II experiment, the total 
accuracy of the results is statistics-limited, with system- 



TABLE III. Total systematic uncertainty in R and its contri- 
butions. See text for details. 



Q^ GeV^ 


3.5 




4.0 




4.8 




5.6 




•&FPP 


1.4 X 10" 




0.8 X 10" 




1.4 X 10" 


r3 — 


0.7 X 10" 


.3 


(fFPP 


5.1 X 10" 


-3 


6.3 X 10" 


-3 


6.1 X 10" 


-3 


2.9 X 10" 


■3 


^bend 


4.6 X 10" 


-3 


0.1 X 10" 


-3 


2.6 X 10" 


-3 


4.3 X 10" 


■3 


't>be.nd 


1.3 X 10" 


-3 


1.1 X 10" 


-3 


6.1 X 10" 


-3 


12.3 X 10^ 


-3 


COSY 


0.4 X 10" 


-3 


0.4 X 10" 


-3 


1.2 X 10" 


-3 


12.7 X 10" 


-3 


Bckgr. 


N. A. 




0.9 X 10" 


-3 


1.1 X 10" 


-3 


1.2 X 10" 


■3 


Cuts 


N. A. 




5.4 X 10" 


-3 


7.2 X 10" 


-3 


3.9 X 10" 


-3 


Total 


7.0 X 10" 


-3 


8.5 X 10" 


3 


11.7 X 10 


-3 


18.9 X 10" 


-3 



atic uncertainties at a much lower level. 

As in the original publication no radiative correc- 
tions have been applied to the data presented here. Stan- 
dard, model-independent radiative corrections to R were 
calculated in Q for kinematics very close to those of the 
GEp-II experiment and found to be less than 1% (rel- 
ative) for all four values. Though even 1% relative 
corrections are much smaller than the statistical uncer- 
tainties in the data, the calculations in f6| were performed 
assuming a much wider "inelasticity" cut than that ef- 
fected by the combination of cuts applied in the present 
analysis, such that in reality, the standard radiative cor- 
rections to the GEp-II data are even smaller, which jus- 
tifies neglecting them here. 

IV. RESULTS 
A. Discussion of the Data 

The final results of the GEp-II experiment are reported 
in Table HV] and presented in Figure [121 The values and 
statistical uncertainties of Pf and Pf presented in Ta- 
ble IIVI (and Figure [TT|) are obtained from equation ([6]) . 
Because the analyzing power is calibrated using equation 
([1]), the extracted Pf and Pf values are, by definition, 
equal to those of equations ([1]), which depend only on 
R and kinematic factors, regardless of the value of Pe 
assumed in the analysis. For reference, the values of Pe 
used in the analysis at each are shown in Table HI 
These values are based on the average of all beam po- 
larization measurements at a given setting. Because Pe 
was stable at the few percent level throughout the du- 
ration of each kinematic setting, a single Pe value was 
assigned to all data taken at a given Q^. As presented, 
the statistical uncertainties in Pf and Pf correspond to 
the uncertainties in the raw asymmetries measured by 
the FPP, which are large compared to the corresponding 
systematic uncertainties. 

P( and Pi can also be calculated from R and kine- 
matic factors using equation ([T]). Neglecting the very 
small covariance of Pj and Pi and the uncertainty in the 
kinematic factors involved, the uncertainty in R is given 
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TABLE IV. Final results of the GEp-II experiment. (Q^) is the acceptance-averaged , while AQ^ is the half-width of the 
total interval, which is centered at the nominal . The raw {Pi''"), background (P?"'^') and corrected (Pf') polarization 
transfer components and the raw form factor ratio 7? are presented with statistical uncertainties only. The background fraction 
(/) averaged over the final cut region is quoted with its systematic uncertainty A/. The final results for R = /ipG^/GJ^ 
are quoted with statistical and systematic uncertainties. The data at = 3.5 GeV^ were not reanalyzed, and the quoted 
result is identical to that of the original publication ^. The originally published results 0] are quoted on the bottom line for 
comparison. 



Nommal Q (GeV ) 


3.5 


4.0 


4.8 


5.6 


(Q")± AQ^ (GeV^) 


3.50 ±0.23 


3.98 ±0.26 


4.76 ± 0.30 


5.56 ±0.34 




N. A. 


-0.108 ±0.011 


-0.094 ±0.011 


-0.070 ±0.017 


pt" ± Ap;*"" 


N. A. 


0.683 ±0.012 


0.795 ±0.013 


0.886 ± 0.030 


R±AR (raw) 


N. A. 


0.514 ± 0.055 


0.445 ± 0.052 


0.350 ± 0.085 


if) ± A/ 


N. A. 


(0.30 ± 0.04)% 


(0.38 ± 0.06)% 


(0.47 ± 0.07)% 


pinel _j_ spinel 


N. A. 


0.116 ±0.051 


0.264 ± 0.038 


0.128 ±0.034 


pinel _j_ pinel 


N. A. 


0.224 ± 0.053 


0.006 ± 0.049 


0.278 ± 0.072 


Pf ± APf 


-0.118 ±0.015 


-0.109 ±0.011 


-0.096 ±0.011 


-0.071 ±0.017 


Pf ± APf 


0.616 ±0.017 


0.685 ±0.012 


0.799 ±0.013 


0.890 ± 0.030 


R ± ARstat ± ARsyst (final) 


0.571 ±0.072 ±0.007 


0.517 ±0.055 ±0.009 


0.450 ± 0.052 ± 0.012 


0.354 ±0.085 ±0.019 


Pt ± APt (equation ([T|)_^ 


-0.118 ±0.014 


-0.109 ±0.011 


-0.096 ±0.011 


-0.071 ±0.017 


Pe ± APe (equation ([T]))3 


0.616 ±0.005 


0.685 ± 0.003 


0.799 ± 0.002 


0.890 ± 0.002 


P ± AP,t,t ± AP,^,i ([2]) 


0.571 ±0.072 ±0.007 


0.482 ±0.052 ±0.008 


0.382 ±0.053 ±0.011 


0.273 ± 0.087 ± 0.028 



These are the values of Ft and Pi calculated from equation JTJ, with uncertainties due solely to the uncertainty in R. 



by {AR./Rf = {APt/Ptf + {/APi/Pif. While the un- 
certainties in Pt and Pg obtained from equation ([6]) are 
similar, and /S.Pf is generally larger than /S.Pf due to 
the unfavorable precession angle, the uncertainty in R is 
nevertheless dominated by the uncertainty in Pj, since 
is generally large compared to Pt- Due to the weak 
sensitivity of to i?, the uncertainty in P^ calculated 
from equation |T]) is much smaller than the uncertainty 
in Pi extracted from the FPP asymmetry. On the other 
hand, since Pt is proportional to R and the relative un- 
certainties in Pt and R are similar, the uncertainty in Pt 
calculated from equation ([1]) is very similar to the uncer- 
tainty in its extraction from the measured asymmetry. 

Figure [12] shows the final results with the GEp-I 
data the orig inally published GEp-II data 0, the 
GEp-III data jig, and the combined GEp-27 result [H], 
consisting of a weighted average of measurements of R 
at three e values for a fixed of 2.5 GeV^. The curves 
illustrate the effect of the revised data on a global fit 
using the Kelly parametrization [42l| of and G\,j to 
elastic ep cross section and polarization data, including 



the GEp-III data 
are detailed in [43 1 



25{ . The data selection and fit method 
The dashed "Old fit" curve uses the 
original GEp-II results, while the solid "New fit" curve 
replaces the three highest-Q^ points from GEp-II with 
the results of the present reanalysis. The combined con- 
tribution of the six highest-Q^ data points to the of 
the "Old" global fit is 2.68. In the "New" fit, the same 
contributior{3 drops to 0.55, indicating a significant 
improvement in the consistency of the data at high . 



The noticeable systematic increase in the results for R 
in the improved data analysis is mostly attributable to 
the systematic effect of the background, underestimated 
by the original analysis. Indeed, the added 5p cut of the 
present work suppresses the background contamination 
to well below 1%, minimizing the associated correction 
and its uncertainty. The most significant difference be- 
tween the final and original results not attributable to 
background or changes in elastic event selection cuts is 
caused by the improved description of the FPP analyz- 
ing power in the present analysis. In the original analysis, 
the momentum dependence of the analyzing power was 
neglected, and the data were divided into discrete bins 
in the FPP polar scattering angle "d. In each bin, the 
analyzing power was extracted from the measured asym- 
metries using equation (|T]) as described in section ITlI B 31 
Then, the analyzing power, which enters equation ([5]) as 
a weight, was assigned to each event in a given ■& bin 
according to the extracted Ay result in that bin. This 
method is approximately equivalent to analyzing the data 
in bins of d assuming Ay = 1, and then combining the 
results of all -d bins in a weighted average to obtain the 
final result. 

In the present work, the final results were obtained 
from a completely unbinned maximum-likelihood analy- 
sis in which Ay{p,pT) was described using the smooth 
parametrization of the pt dependence presented in sec- 
tion [IITB^I and a global momentum scaling Ay cx 1/p, 
leading to a slightly different relative weighting of events 



In this context, the contribution of the six data points is 



defined as J^Ll {^i - RfitiQDY /{^RiY\ 
malized by the number of data points. 
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• JonesOO 




2 4 6 8 

Q2 (GeV^) 

FIG. 12. (color) Polarization transfer data for G^/G^ 
from ;i] (JonesOO), (Gayou02), [H (PuckettlO), ^ 
(Mezianell) and the present work. Error bars are statistical. 
The data of 2] are offset slightly in for clarity. System- 
atic uncertainties for the present work and [25|] are shown as 
bands below the data. The inset shows a zoomed view of 
the data near — 2.5 GeV^, demonstrating the excellent 
agreement between high-precision data from Hall A [3| and 
Hall C ^26] at this Q^. Curves are global proton form factor 
fits using the orig inally published GEp-II data 0] (Old fit) 
and the present work (New fit), with standard la point- wise 
uncertainty bands. Both fits include the GEp-III data. The 
linear fit of equation ((2| is shown for comparison. See text 
for details. 



as a function of p and d. At = 5.6 GcV^, where 
the statistical uncertainty is large, roughly half the dif- 
ference between the originally published result and the 
final result is attributable to the different description of 
Ay (with the other half coming from the background), 
while at = 4,0 and 4.8 GeV^, the effect of the Ay 
description is small and the difference is dominated by 
the background effects. This observation can be under- 
stood by examining the 1) dependence of R in Figure [5] 
and the pr-dependence of Ay in Figure IH] at = 5.6 
GeV^. A negative fluctuation of R in the d bin near 11° 
coincides with a positive fluctuation of Ay in the pT bin 
near 0.7 GeV. Assigning this value of Ay to all events in 
this bin artificially overweights the corresponding nega- 
tive statistical fluctuation in R, inducing a slight negative 
bias to the result. Because this particular fluctuation is 
relatively large, the effect of using a smooth parametriza- 
tion of the analyzing power instead of a discretely binned 
description is noticeable. This is in contrast to the two 
lower-Q^ points, for which no large t?-dependent statis- 
tical fluctuations of Ay or R are observed, making the 
combined result rather insensitive to the description of 



Ay. It cannot be too strongly emphasized that the de- 
pendence of the result on the description of Ay derives 
only from p and i?-dependent statistical fluctuations of 
the data, since Ay cancels in the ratio Pt/Pi (see Fig- 
ure [8] and the discussion in section IIII B 3p . Therefore, 
the sensitivity of the results to the description of Ay is 
properly regarded as part of the statistical uncertainty, 
and no additional systematic uncertainty contribution is 
assigned. 

Despite discarding up to 10% of the events included in 
the original analysis, the statistical error of the final re- 
sult for R is actually slightly reduced at = 4.8 and 5.6 
GeV^ relative to the original publication. The improve- 
ment reflects an increase in the effective Ay of the final 
sample of events due to the improved suppression of the 
background, which tends to dilute the measured asymme- 
try. On the other hand, the statistical error at = 4.0 
GeV^ has slightly increased relative to the original pub- 
lication, since at this the loss of statistics slightly 
outweighs the increase in Ay from improved background 
suppression. Nonetheless, the quality of the result is im- 
proved by the removal of a previously underestimated 
systematic error. 

Compared to the situation before the GEp-III exper- 
iment, the emerging picture of the large-Q^ behavior of 
G^^/G\,j is considerably modified. Before GEp-III, the 
GEp-I and GEp-II data suggested a strong linear de- 
crease of R continuing to high . The linear trend of 
the data suggested a zero crossing of G^^/G\j before 8 
GeV^. The GEp-III data showed that the linear decrease 
probably does not continue to higher Q^, at least not at 
the slope suggested by the GEp-I and original GEp-II 
results. Although the lower-Q^ data from GEp-27 ap- 
peared to rule out any neglected systematic error in the 
GEp-III data, the fact that all three data points from 
GEp-III were systematically above the trendline of the 
previous data raised concern about the consistency be- 
tween different experiments and the reproducibility of 
the polarization transfer method. Moreover, while there 
was no a priori reason to expect the linear decrease to 
continue, and the apparent ~ l.Scr disagreement between 
GEp-II and GEp-III did not rise to the level of statistical 
signiflcance, the lessons learned from the GEp-III analy- 
sis, particularly with respect to backgrounds, motivated 
a reanalysis of the GEp-II data, leading to the results 
presented in this article. With imp roved analysis, the 
data from Halls A and C [l|, H, l25l [26| are now in ex- 
cellent agreement over a wide range, bringing added 
clarity to the experimental situation regarding G^/G^j. 

In a simple global analysis using the Kelly parametriza- 
tion [4^, the data before GEp-III implied a zero cross- 
ing at = 9 GeV^, with an uncertainty range of 
7.7 GeV^ < < 12.5 GeV^, based on the pointwise 
Icr error bands of the fit result. After adding the GEp- 
III data and replacing the GEp-II data with the final 
analysis results, the zero crossing is shifted to 15 GeV^, 
with an uncertainty range of roughly 12 GeV^ < < 29 
GeV'^. Although the size of the error band in G^ shrinks 
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by a factor of two at large when the GEp-III data 
are added, the reduced slope of G\ increases the uncer- 
tainty in the location of the potential zero crossing. The 
Kelly parametrization, despite having the correct static 
limit and sensible pQCD-based asymptotic behavior at 
high Q^, does not describe the actual physics involved 
in the transition between low and high-Q^ asymptotic 
behavior. Therefore, its extrapolation beyond the range 
of the existing data necessarily understates the true un- 
certainty in the behavior of at large Q^. Only future 
measurements at higher with higher precision can 
definitively reveal the behavior of in the region where 
the predictions of leading models of the nucleon diverge, 
as discussed in the following section. 



B. Physics Interpretation 

1. Perturbative QCD 

Perturbative QCD (pQCD) makes rigorous predictions 
for the dependence of the nucleon form factors when 
is sufficiently large that the scattering amplitude can 
be factorized as the convolution of a baryon distribution 
amplitude with a perturbatively calculable hard scatter- 
ing kernel [i^. At leading order in the Dirac form 
factor is proportional to as/Q"^ times slowly varying loga- 
rithmic terms, because the large momentum transfer ab- 
sorbed by the struck quark must be shared among the two 
spectator quarks via two hard gluon exchanges in order 
for the nucleon to recoil as a whole. The Pauli form fac- 
tor is suppressed by an extra power of at leading order 
due to helicity conservation implying that Q^F2/Fi 
(and therefore Ge/Gm) becomes constant at very high 

. While pQCD predicts the asymptotic dependence 
of the nucleon form factors, it does not predict the value 
of at which the hard scattering mechanism becomes 
dominant. Isgur and Llewellyn Smith [46l l47j have ar- 
gued that pQCD is not applicable to observables of exclu- 
sive reactions such as the nucleon form factors in the ex- 
perimentally accessible region. Ralston and Jain [i^ , 
inspired by the results of the GEp-I and GEp-II exper- 
iments, revisited the leading power behavior in l/Q of 
F2/F1 in the pQCD hard-scattering picture by consider- 
ing the violation of hadron helicity conservation that en- 
sues when quark wavefunction components with nonzero 
orbital angular momentum are included, and found that 
F2/F1 cx l/Q. 

Belitsky, Ji and Yuan like Ralston and Jain '48'], 
argued that quark orbital angular momentum is the dom- 
inant mechanism for nucleon helicity flip at large 
in pQCD, owing to the very small mass of the current 
quarks involved in the hard scattering. They performed 
a pQCD analysis of the proton's Pauli form factor F2 in- 
cluding the subleading-twist contribution to the proton's 
light-cone wavefunction. The leading-order pQCD con- 
tribution to F2 involves intial and final-state light-cone 
wavefunctions differing by one unit of quark orbital an- 
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FIG. 13. (color online) Data for the ratio /Ff from selected 
polarization experiments including Refs. 1, 25, 32, 49-53] and 
the final GEp-II data from Table IIVI of this work. From top 
to bottom: Q^F^/F[, QF^/Ff and QVln^ (Q'^/J^^) F^/Ff, 
where A = 236 MeV was fitted to the data for > I GeV^ 
The curves in the top panel are GPD model fits from [13] 
(Guidal05) and [3 (Diehl05). 



gular momentum, with zero orbital angular momentum 
in either the initial or final state. In this calculation, log- 
arithmic singularities in the convolution integrals lead to 
the modified scaling Q^F^/F^ cx In^ (Q'^/A^), where A 
is an infrared cutoff parameter related to the size of the 
nucleon. 

Figure [13] shows the experimental data for F2 /Ff plot- 
ted as Q^Fi/Ff, QFi/Ff and Q^/W (QVA^) Ff/f'f. 
Clearly, the leading- twist, leading-order pQCD scaling 
behavior is not respected by the data in the presently 
accessible region, although the slope of 
does appear to be trending toward a fiat behavior at 
the highest-Q^ values measured so far. The scaling of 
QF2/F1 predicted by is approximately satisfied up 
to 8.5 GeV^, although there is a hint that F2 may start 
to fall faster than Fi/Q for higher Q. The logarithmic 
scaling of [12] is satisfied for ^ 1 GeV^ at a value of 
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the cutoff parameter A = 236 MeV [hc/k = 0.835 fm) 
determined by fitting the data for > 1 GeV^. 

While the "precocious" scaling of F2 / Ff is interesting, 
it is probably largely accidental, perhaps a consequence 
of delicate cancellations of higher-order effects in the ra- 
tio [l^l ■ The scaling of F2 /Ff is a necessary but insuffi- 
cient condition for the onset of the perturbative regime. 
pQCD-based form factor predictions based on light-cone 
sum rules [13, [13 have yet to reach the level of accu- 
racy achieved by the phenomenological models discussed 
below in describing all four nucleon form factors. In the 
GPD model fits shown in Figure[T31 the 'Feynman' mech- 
anism corresponding to the overlap of soft wavefunctions 
dominates the form factor behavior. The neutron data 
for F2 /F" do not scale in the currently measured re- 
gion up to 3.4 GeV^ for values of A similar to that which 
describes the proton data [2^. Moreover, combining the 
proton and neutron data to separate the up and down- 
quark contributions to the nucleon form factors ;2^] re- 
veals that the ratios f^/F" and F^/F^ become approx- 
imately constant above 1 GeV^, at odds with the asymp- 
totic pQCD picture, while the ratios Ff / F^ and F^ / F^ 
decrease at high Q^, a behavior that can be explained in 
terms of diquark degrees of freedom [T^ . Based on these 
and other considerations, it is generally believed that the 
nucleon form factors are dominated by non-perturbative 
physics in the 1-10 GeV^ region addressed by present 
experiments. 



2. Generalized Parton Distributions 

The generalized parton distributions (GPDs) are uni- 
versal non-perturbative matrix elements involved in the 
QCD factorization of hard exclusive processes such as 
deeply virtual Compton scattering (DVCS) [ll[li-[li]. 
The GPDs are functions of the longitudinal momentum 
fraction x, the momentum fraction asymmetry or "skew- 
ness" ^ and the squared momentum transfer to the nu- 
cleon t (not to be confused with the photon virtuality 
Q^). GPDs play a crucial role in the synthesis of seem- 
ingly disparate nucleon structure information obtained 
from inclusive and exclusive reactions. The Dirac and 
Pauli form factors Fi and F2 equal the first x moments 
of the vector {H{x,t)) and tensor {E{x,t)) GPDs, re- 
spectively. In the forward {t — 0) limit, H{x,t = 0) is 
the valence quark density. Precise measurements of the 
Pauli form factor F2 at large constrain the behavior 
of E(x^ t), yielding new information on nucleon structure 
that is inaccessible in inclusive deep inelastic scattering 
(DIS). With increasing Q^, the strength in the GPD in- 
tegrals corresponding to the form factors is increasingly 
concentrated in the high- a: region. Therefore, the a; — > 1 
behavior of H{x,t) and E{x,t) can be constrained by 
fitting the high-Q^ nucleon form factors. 

While systematic studies of the observables of DVCS 
and other hard exclusive reactions promise an eventual 
direct extraction of GPDs from global analysis (for re- 



cent examples, see [59|-|63j), the experimental mapping 
of these observables is still at an early stage. Meanwhile, 
constraints from the elastic form factors and the forward 
parton distributions measured in DIS have been explored 
using physically motivated GP D p arametrizations based 
on Regge phenomenology [3, [ll|. In both models, the 
high- a: behavior of E was determined by the high 
behavior of i^f measured by the GEp-I and GEp-II ex- 
periments, enabling an evaluation of Ji's sum rule [l3l[57| 
for the total angular momentum carried by the up (J") 
and down {J'^) quarks in the nucleon. The calculations 
of [13 found 2J'^ = -0.06 and 2 J" = +0.58, in qualita- 
tive agreement with lattice QCD calculations available at 
the time 64 1 , as well as more recent calculations [6M [l^ . 
The predictions of the GPD models of [T2| and fl5| are 
compared to the data for Q^F2 /Ff in Figure [T51 

The two-dimensional Fourier transform of the t- 
dependence of the GPDs at ^ = yields a three- 
dimensional impact parameter representation p{x,h±) 
in two transverse spatial dimensions and one longitudi- 
nal momentum dimension [67| . By forming the charge- 
squared weighted sum over quark fiavors and integrating 
over all x, Miller [l6| derived the model-independent in- 
finite momentum frame transverse charge density Pch{b) 
as the two-dimensional Fourier transform of the Dirac 
form factor Fi. The Pauli form factor F2 can also be 
related to the transverse anomalous magnetic moment 
density /9m(b) [11]. Miller et al. [1^ performed the first 
analysis of the uncertainties in the transverse charge and 
magnetization densities of the proton due to the uncer- 
tainties and incomplete coverage of the form factor 
data. Measurements of at yet higher are needed 
to reduce the uncertainty in Pm{b) at small b. 

Since an exact solution to QCD in the non- 
perturbative regime is not yet possible, predicting nu- 
cleon form factors in the domain of strong coupling and 
confinement is rather difficult. Consequently, many phe- 
nomenological models have aimed to unravel the compli- 
cated internal structure of the nucleon in this domain. 
The following discussion provides an overview of a wide 
range of models. 



3. Vector Meson Dominance 

The global features of the nucleon form factors were 
explained by early models based on vector meson dom- 
inance (VMD) 170|. VMD models are a special case 
of dispersion relation analyses, which provide a model- 
independent, non-perturbative framework to interpret 
the electromagnetic structure of the nucleon in both the 
space-like and time-like regions. Early VMD model cal- 
culations included the p and its excited states for the 
isovector form factors, and the u) and for the isoscalar 
form factors. The number of mesons included and the 
coupling constants and masses can be varied to fit the 
data. In practice, many parameters are fixed or strongly 
constrained by experimental data, including but not lim- 
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ited to nucleon FF data, reducing the number of free 
parameters and increasing the predictive power of the 
approach. More recent calculations have used the pQCD 
scaling relations to constrain the large behavior of the 
fits. An example is Lomon's fit |71|], which uses p(770), 
a;(782), ?!)(1020), p'(1450) and w'(1420) mesons and has 
a total of 12 variable parameters [7l|, [t^]- Bijker and 
lachello 73] updated the 1973 model of lachello, Jack- 
son and Lande [TOj . performing a new fit including the 
p(770), a;(782), and 0(1020) mesons, and a phenomeno- 
logical "direct coupling" term attributed to an intrinsic 
three-quark structure of rms radius ^ 0.34 fm. 

Despite the relatively good fits obtained by VMD mod- 
els, the approach is at odds with general constraints from 
unitarity. This difficulty can be overcome using disper- 
sion relations. Hohler's dispersion relation analysis [t^] 
was extended in the mid-nineties by Mergell, Meissner, 
and Drechsel [tH] to include nucleon form factor data in 
the time-like region [t^. The analysis of 75 1 has been 
further improved by Belushkin tt al. (77j . In addition 
to the 2tt continuum present in the isovector spectral 
functions, the pn and KK continua were included in the 
isoscalar spectral functions. In [t^], the 2tt continuum 
was reevaluated using the latest experimental data for 
the pion form factors in the time-like region. A simul- 
taneous fit to the world data for all four form factors in 
both the space-like and time-like regions was performed. 
The results are in very good agreement with the data 
available at the time. Dubnicka et al. developed a uni- 
tary and analytic ten-resonance model including the 27r 
continuum [tsI, [tII , which fits all nucleon form factors in 
both the space-like and time-like regions. 

Figure [U compares the predictions of selected VMD- 
based models to the experimental data for /ipG^/G^^. 
Of the models shown, the latest version of Lomon's fit 
[T^ I with twelve adjustable parameters achieves the best 
overall agreement with the data for all four form factors 
at spacelike Q^, emphasizing a smooth evolution from 
VMD behavior at low to pQCD scaling at asymptot- 
ically high . Apart from fitting to a more complete 
data set, the main added feature of the model of Bijker 
and lachello [t^ relative to the 1973 model of lachello, 
Jackson and Lande is the inclusion of a "direct" coupling 
term in the isoscalar Pauli form factor which improves the 
large-Q^ behavior of and G^. This model achieves a 
rather good fit to all four FFs using just six adjustable 
parameters (compared to five in the 1973 model). 




5 10 

Q2 (GeV^) 



FIG. 14. (color online) Comparison of selected VMD model 
predictions from ^ (Iachello73), ^ (Lomon06), ^ (Bi- 
jker04), and [t^I (AdamuscinOS) to selected fipG^/G^j data. 
Data are from cross section [s^, [U (empty circles) and 
polarization [l|, (filled circles) measurements, where the 

results of 0] are replaced by the results of the present work. 
See text for details. 



recent calculations focus on the isovector form factors, 
for which the contributions from disconnected diagrams 
are reduced. Calculations are performed for various TOtt 
values and lattice spacings a and then extrapolated to 
the physical pion mass and the continuum limit a ^ 0. 
Recently the QCDSF/UKQCD Collaboration has per- 
formed calculations j8^] at = 180 MeV with different 
lattice spacings and volume sizes, but the upper range 
is limited to 3 GeV^. Lattice QCD form factor calcula- 
tions in the region measured by the GEp-II and GEp- 
III experiments are difficult due to large statistical and 
systematic errors. Calculations by Lin et al. employ a 
novel technique to extend the reliable range of the cal- 
culations to = 6 GeV^ at m,r > 450 MeV for quenched 
and dynamical ensembles [s^. Nonetheless, calculations 
at such high must ultimately be performed with a 
finer lattice spacing to reduce the systematic error. 



4. Lattice QCD 

Lattice QCD calculations provide ab initio evaluations 
of static and dynamic hadron properties, including the 
nucleon electromagnetic form factors, from numerical so- 
lutions of QCD on a finite-volume lattice of discrete 
space-time points. At present, the lattice calculations 
are done using unphysically large quark masses which 
are given in terms of the pion mass, ttitt. Moreover, most 



5. Constituent Quark Models 

In the constituent quark model (CQM), the nu- 
cleon consists of three constituent quarks, which can 
be thought of as valence quarks that become much 
heavier than the elementary quarks appearing in the 
QCD Lagrangian when dressed by gluons and quark- 
antiquark pairs. The dressing effects are absorbed into 
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the masses of these quasiparticle effective degrees of free- 
dom. The early non-relativistic constituent quark model 
achieved considerable success in describing the spectrum 
of baryons and mesons with correct masses [13] . In or- 
der to describe dynamical quantities such as form factors 
in terms of constituent quarks, a relativistic description 
(rCQM) is mandatory because the values involved in 
modern experiments have reached as high as ten times 
the nucleon mass squared and ~ 10^ times the "bare" 
quark mass squared. 

Frank et al. |85| calculated and G^j in the light- 
front CQM using the light-front nucleon wave function of 
Schlumpf [s^, and predicted that G^ might change sign 
near 5.6 GeV^, a behavior inconsistent with current data, 
though qualitatively correct. In this model, constructing 
a Poincare-invariant nucleon wavefunction that is also 
an eigenstate of spin leads to the substantial violation 
of hadron helicity conservation [87] responsible for the 
observed scaling of QF2/F1 in the range of present 
experiments. This feature is a consequence of the unitary 
Melosh rotation 'ss'] which mixes quark spin states in the 
process of boosting the nucleon spin-flavor wavefunction 
from the rest frame to the light front. Miller extended 
this model to include pion-cloud effects important 
to the understanding of the low-Q^ behavior generally 
and G^ in particular. 

Gross et al. [q^, IU modeled the nucleon as a bound 
state of three dressed valence constituent quarks in the 
covariant spectator formalism, in which the virtual pho- 
ton is absorbed by an off-shell constituent quark, and 
the two spectator quarks always propagate as an on-shell 
diquark. In this model, the constituent quarks have in- 
ternal structure described by form factors which become 
pointlike at large as required by pQCD and exhibit 
VMD-like behavior at low . The model nucleon wave- 
function of ^91] obeys the Dirac equation, includes only s- 
wave components, and its spin-isospin structure reduces 
to that of the SU(2) x SU(2) quark model in the non- 
relativistic limit. 

Cardarelli et al. calculated the ratio using light- 
front dynamics and investigated the effects of SU(6) sym- 
metry breaking. As in [s^, they showed that the decrease 
of R with increasing is caused by the relativistic effect 
of the Melosh rotations of the constituent quark spins. 
De Sanctis et al. calculated the nucleon form factors 
in the relativistic hypercentral constituent quark model 
(hCQM) [9^. A good fit to all the nucleon form factors 
was obtained using a linear combination of monopole and 
dipole constituent quark form factors. The calculation 
was recently extended to = 12 GeV^ (9J|]. The same 
group also performed calculations within the relativis- 
tic interacting quark-diquark model [o^ , which does not 
achieve the same level of agreement with the data as the 
hCQM. 

De Melo et al. [9^ examined the non-valence compo- 
nents of the nucleon state in light-front dynamics, achiev- 
ing a good description of all spacelike and timelike nu- 
cleon FF data with the inclusion of the Z-diagram involv- 
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FIG. 15. (color online) Selected rCQM predictions from 
Pi (Miller02), Jgl (Gross08), ^ (CardarelliOO), @ (San- 
topintolO) and ^ (deMelo09) for R = Jh,G^/Glj, compared 
to selected data from cross section [^, [sol . [sH (empty circles) 
and polarization P, 0, (filled circles) experiments, wiiere 
the results of 01 are replaced by the results of the present 
work. See text for details. 



ing qq pair creation in addition to the triangle (valence) 
diagram. The chiral constituent quark model based on 
Goldstone-boson-exchange dynamics was used by BofE. et 
al. 97] to describe the elastic electromagnetic and weak 
form factors in a covariant framework using the point- 
form approach to relativistic quantum mechanics. 

Figure [T5l compares the predictions of selected rCQM 
calculations to selected data for R. Of the calcu- 
lations shown, those in which the constituent quarks 
have internal structure represented by CQ form factors 
([m, m, |9i]) and/or significant VMD-related contribu- 
tions to the photon- nucleon vertex (^]) describe the 
data better than those in which the constituent quarks 
are pointlike ([89*]) and have only direct coupling to the 
photon. Although this may be related to the greater 
number of adjustable parameters in models with CQFFs, 
it is apparently physically meaningful that most of the 
models require structure of the constituent quarks and/or 
significant nonvalence {qq pair creation) contributions to 
achieve a good description of the data. 



6. Dyson-Schwinger Equations 

A different theoretical approach to the prediction of 
nucleon form factors is based on QCD's Dyson-Schwinger 
equations (DSEs). The DSEs are an infinite tower of 
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coupled integral equations for QCD's Green functions 
that provide access to emergent phenomena of non- 
perturbative QCD, such as dynamical chiral symme- 
try breaking and confinement 98]. The DSEs admit a 
symmetry-preserving truncation scheme that enables a 
unified description of meson and baryon properties. The 
approach has already achieved considerable success in 
the pseudoscalar meson sector [19|j. The prediction of 
nucleon form factors in the DSE approach involves the 
solution of a Poincare-covariant Faddeev equation. In 
the calculations of dressed quarks form the elemen- 
tary degrees of freedom and correlations between them 
are expressed via scalar and axial vector diquarks. The 
only variable parameters in this approach are the diquark 
masses, fixed to reproduce the nucleon and A masses, 
and a diquark charge radius embodying the electro- 
magnetic structure of the diquark correlations. A dif- 
ferent approach to DSE-based form factor calculations 
effects binding of the nucleon through a single dressed 
gluon exchange between any two quarks (isj without ex- 
plicit diquark degrees of freedom. In this calculation, the 
only parameters are a scale fixed to reproduce the pion 
decay constant and a dimensionless width parameter i] 
describing the infrared behavior of the effective coupling 
strength of the quark-quark interaction. 



The predictions of several DSE-based calculations for 
the proton Sachs form factor ratio R = f^pG^^/G^j^j are 
shown in Figure 1161 The quark-diquark model calcula- 
tion [l3| underpredicts the data at low but agrees rea- 
sonably well at higher . The disagreement at low is 
attributed to the omission of meson cloud effects. The ad- 
dition of dynamically generated, momentum-dependent 
dressed-quark anomalous magnetic moments that be- 
come large at infrared momenta improves the description 
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of R at low Q^. The three-quark model calculation 
agrees with the data at low Q^, but underpredicts the 
data at higher Q^, becoming numerically unreliable for 
Q2 > 7 GeV2. 



The deficiencies of the DSE approach, including the ap- 
proximation schemes required to make the calculations 
analytically tractable and the omission of meson-cloud 
effects, are evident in the disagreement between the pre- 
dicted form factors and the experimental data, which is 
more severe than in the various models described above, 
which have more adjustable parameters. The advantage 
of the approach is that it provides a systematically im- 
provable framework for the ab initio evaluation of hadron 
properties in continuum non-perturbative QCD, that is 
complementary to discretized lattice simulations. As fun- 
damental measurable properties of nucleon structure, the 
electromagnetic form factors are essential to the feedback 
between theory and experiment required to make further 
progress in this direction. 
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FIG. 16. (color online) Predictions of DSE-based calcula- 
tions for R = lJ.pG^/G^ compared to experimental data 
from cross section [a, ISCl . IsH (empty circles) and polariza- 
tion [H, 0, [2^ (filled circles) experiments, where the results of 
^ are replaced by those of the present work. The results of 
13 (CloetOQ) are shown for a particular choice of the diquark 
charge radius. The curve from [O^ (Changll) is that of [13] 
with the addition of dressed quark anomalous magnetic mo- 
ments. The results of ^1§\ (Eichmannll) are shown for two 
values of 77, showing the weak sensitivity of the form factor 
results to this parameter. 



7. AdS/QCD 



In the past decade, theoretical activity has flourished 
in modeling QCD from the conjecture of the anti-de 
Sitter spac e / conform al field theory (AdS/CFT) corre- 
spondence 133l - ll35 |. a mapping between weakly coupled 
gravitational theories in curved five-dimensional space- 
time and strongly coupled gauge theories in flat four- 
dimensional space-time. Since QCD is not a conformal 
field theory, the symmetry of the anti-de Sitter space 
is broken by applying a boundary condition. Brodsky 
and de Teramond 136] have calculated Fi for the pro- 
ton and neutron and emphasized the agreement of the 
predicte d Q^ Fi dependence with the data. Abidin and 
Carlson [137l | have calculated both proton and neutron Fi 
and F2 along with the tensor form factors using both hard 
and soft wall boundary conditions. This model predicts 
the same asymptotic dependence as the dimensional 
scaling of pQCD, but does not reproduce the detailed 
features of the data in the presently measured region. 
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FIG. 17. (color online) Comparison of selected theoretical predictions to data for all four nucleon FFs at space-like Q^. Theory 
curves are 15] (DiehlOS), 18] (Eichmannll), ^ (Lomon06), [ol (GrossOS) and [93 | (SantopintolO). data are from 
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8. World nucleon form factor data compared to theory 

Figure [17] summarizes the theoretical interpretation 
of the nucleon electromagnetic form factors, with rep- 
resentative examples from each of the classes of models 
discussed compared to the world data for all four nu- 
cleon electromagnetic form factors. Published results for 
R — iipG^^/G\j were converted to values using the 
global fit of and G\,j from [isj , updated to use the R 
values of the present work, a change that does not notice- 
ably affect G\j. Except at very low , the contribution 
of the uncertainty in G^^ to the resulting uncertainty in 
G^ is negligible. At this juncture, it is worth recalling 
that the G^ results extracted from cross section data are 
believed to be unreliable at high Q"^ due to incompletely 
understood TPEX corrections, which have not been ap- 
plied to the data shown in Figures [HlfTTl Except for the 
DSE calculation of all of the models shown describe 



existing data very well, which is to be expected given that 
the parameters of the models are fitted to reproduce the 
data. However, their predictions tend to diverge when 
extrapolated outside the range of the data. That 
the DSE-based calculation of fails to describe the 
data as well as the other calculations is not surprising, 
since it represents a more fundamental ah initio approach 
with virtually no adjustable parameters, but requires ap- 
proximations that are not yet well-controlled. Significant 
progress in the quality of the predictions is nonetheless 
evident, as the data expose the weaknesses of different 
approximation schemes. Since the hard scattering mech- 
anism leading to the asymptotic pQCD scaling relations 
is not expected to dominate the form factor behavior at 
presently accessible values, phenomenological mod- 
els and the ambitious ongoing efforts in lattice QCD and 
DSE calculations are of paramount importance to under- 
standing the internal structure and dynamics of the nu- 
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cleon. Planned measurements at higher following the 
12 GeV upgrade of JLab promise to be of continuing in- 
terest and relevance owing to their power to discriminate 
among the various models and to guide the improvement 
of the more fundamental calculational approaches. 

V. CONCLUSION 

This work has presented an expanded description and 
an improved final data analysis of the GEp-II experiment, 
originally published in [2], which measured the proton 
electromagnetic form factor ratio for 3.5 GeV^ < < 
5.6 GeV^ in Jefferson Lab's Hall A using the polariza- 
tion transfer method. The improved data analysis finds 
a systematic increase in the results for R = /ipG^/G^^ 
that improves the agreement between the GEp-II and 
GEp-III 25] data. This increase mainly reflects the un- 
derestimated impact of the tt*^ production background in 
the original analysis of GEp-II. Section |ll] presented the 
details of the experimental apparatus and described the 
differences between the GEp-II and GEp-III experiments. 
Section [Till presented the full details of the data analysis, 
including the selection of elastic events in section IIII Al 
the extraction of polarization observables in section lTlIBl 
and the treatment of the background in section lTlI CI The 
analysis of systematic uncertainties was presented in sec- 
tion IIIIDI In section IIV Al the features of the data and 
the sources of the increase in the results relative to the 



original analysis were discussed at length. An overview 
of recent progress in the theoretical understanding of nu- 
cleon form factors was given in section IIVBI In conclu- 
sion, this work represents the final results of the GEp-II 
experiment. The revised data presented here and the re- 
sults of the GEp-III experiment 25] have considerably 
improved the experimental knowledge of at large . 
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